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ABSTRACT 

The starting transients of high-performance solid-pro- 
pellant rocket motors having small port-to-throat area ratios 
are inaccurately predicted and unsatisfactorily analyzed by 
either the widely used lumped chamber -parameter models, or 
the quasi-steady one-dimensional flow models. In this study, 
a comprehensive analytical model has been developed to describe 
the transfer of mass, momentum and heat, and the fast ignition 
propagation along the port of monolithic rocket motors with 
high internal gas velocities during the three major phases 
of the starting transient: induction, flame spreading, and 

chamber filling. Both time and space development of the flow 
field in the motor are considered. 

The gas dynamics in the chamber is governed by a set of 
three hyperbolic partial differential equations, that are 
coupled with the ignition and flame spreading events, and 
with the axial variation of mass addition. The flame spread- 
ing rate is calculated by local successive heating-to-igni- 
tion along the propellant grain surface. Local augmentation 
of the burning rate due to erosive burning is considered by 
using a modified Lenoir-Robillard 1 s law. Entrance, bound- 
ary-layer, and surface roughness effects are accounted for 
in the empirical expressions for heat-transfer and friction 
coefficients . 

The numerical solution of the governing equations 
utilizes an implicit finite-difference scheme combined with 
suitable initial conditions and proper physical and extran- 
eous boundary conditions at both ends of the motor. The 
solution yields pressure, velocity and temperature distri- 
butions along the port for each calculation time step. 

Experimental diagnostic studies were carried out with 
a laboratory-size rectangular window motor (50 cm grain 
length, 5 cm burning perimeter, and 1 cm hydraulic port 
diameter) , using a controllable, head-end gaseous pyrogen 
igniter. Pressure measurements were taken at five stations 
along the port. Tests were conducted with unmetallized AP 
composite propellant at port-to-throat area ratios (Ap/A t ) 
of 2.0, 1.5, 1.2, and 1.06. Accordingly, the head-ena 
pressures varied from 37 to 75 atm. 

Calculated pressure transients tp(t,x)] are in very 
good agreement with those measured in the experimental 
system, for all port-to-throat area ratios tested. Agree- 
ment within 10% has been obtained between calculated and 
experimentally measured ignition-delay and flame-spreading 
times, pressure peaks and times to achieve the pressure peaks 
for Ap/At = 2.0, 1.5 and 1.2. For Ap/At = 1*06, the agree- 
ment in times was within 15%, due to high gas velocities and 
low induction pressures. The highest pressurization rahes 



in all tests were measured at mid-motor. Rates as high as 
7,000 atm/'sec with Ap/A^ = 2.0 and as low as 3,000 atm/sec 
with Ap/A t = 1.06 were measured. 

Parametric studies have shown that the induction period 
and flame spreading duration are greatly affected by the ig- 
niter mass flow rate and gas temperature, the correlation for 
heat-transfer coefficient, Ap/A t between 1.06 and 1.5, and 
uncertainty in the thermal conductivity and ignition temper- 
ature of the propel. ant. The maximum chamber pressure and 
pressurization rate are affected by the igniter mass flow 
rate, burning surf ace-to-throat area ratio, burning rate law, 
Ap/A t , and axial distance along the port. 

The analysis enables starting thrust transient predic- 
tion and control by proper calculation of the stagnation 
pressure at the nozzle entrance. 

An important area for future work is to extend this 
analysis to segmented motors. 

The formulation of the analytical model allows an easy 
extension of the analysis to various types of solid propel- 
lants and operational motor configurations. 
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NOMENCLATURE 

2 

A^ = burning surface area of propellant, cm 

2 

A = cross-sectional area of motor, cm 
m 

2 

Ap = cross-sectional area of the port (port area) , cm 

2 

A fc = motor nozzle throat area, cm 
A^,A 2 ,A 3 = points in Figs. 8 and 9 

a = pre-exponential factor in the non-erosive burning 
rate law (r Q = ap n ) 

a. = coefficient in the expression for convective heat- 
1 transfer coefficient [see Eq. (III-27)] 

a^,a 2 ,a 3 ,a i j,a 5 ,a, = coefficients in correlations and cor- 
x 3 rection functions for the convective 

heat-transfer coefficient (see Section 
D, Chapter III) 

B S/ B r ~ points at left boundary and right boundary, res- 
pectively (see Figs. 8 and 9) 

b = burning perimeter, cm 

c = speed of sound, cm/sec; when with subscript - 
specific heat, cal/g-°K 

c* = characteristic velocity, cm/ sec 

c p = specific heat at constant pressure, cal/g-°K 

d = duct diameter, cm 

d h = hydraulic diameter of the port, cm 

dp = diameter of a cylindrical port, cm 

2 

E = total stored energy per unit mass, e + u /2g«^, cal/g 

E = activation energy of propellant surface reaction, 
cal/g-mole 

e = internal energy of gas, cal/g 

F = functional coefficient [see Eq. (IV-2)] 

2 

f = friction coefficient, 2gt w /pu 


IX 


= volumetric loading-density factor [see Eq. (1-3)] 

4 

= mass velocity, pu/A , g/cm -sec 

ir 

- acceleration of gravity, conversion factor, 
g-cm/gf-sec * 

= lo^al convective heat-transfer coefficient, cal/ 
cm -sec-°K 

= local convective heat-transfer coefficient over the 
propellant surface, cal/cm -sec-°K 

= local convective heat-transfer coefficient over the 
non-propellant port wall, cal/cm -sec-°K 

= enthalpy of the propellant combustion gas at the 
adiabatic flame temperature, cal/g 

= port height (see Fig. 2) , cm 

= specific impulse of the solid propellant, gf-sec/g 

= total impulse delivered by a rocket motor, kgf-sec 


I u ,i T ,I p 


inhomogeneous terms in the governing equations 
[see Eq. (IV-2)] 


J = mecnanical equivalent of heat, gf-cm/cal 
c 

k = erosive burning coefficient, cm J -°K/cal 
k^k^k^k^ = coefficients in Eq. (III-40) 

L = effective length of motor combustion chamber, cm 
L* = characteristic length, cm 

M = Mach number 

m * mess velocity exponent in correlations for the con- 
vective heat-transfer coefficient (see Section A, 
Chapter II) 

m ig * i9 n ^ tei * mass flow rate, g/sec 

m n = mass flow rate through the motor nozzle, g/sec 

N * number of spacewise mesh points for the numerical 
calculation 

Nu « Nusselt number 


*gf denotes unit of gram force. 




X 


n = pressure exponent in the non-erosive burning rate 
law (r Q * ap n ) 

Pr = Prandtl number 

P w = wetted perimeter of the port, cm 

2 

p = static pressure, gf/cm 

2 

p_ = ambient pressure, gf/cm 

aL 

p = stagnation pressure at the motor aft end (nozzle 
entrance) , gf/cnr 

2 

Pind = pressure level during the induction interval, gf/cm 

2 

P t = static pressure at the motor nozzle throat, gf/cm 

p ir p 2 ,p 3 ,p 4 ,p_ = static pressure at stations 1 through 5, 

respectively (see Fig. 2) , gf/cin 

. 2 
p^ = pressure at station 1 at time t^ g , gf/cm 

2 

= rate of heat transfer, cal/cm -sec 

R = specific gas constant of the combustion gases, 
gf-cm/g-°K 

Re = Reynolds number 

R u * universal gas constant, 1.98 cal/g-mole-°K 

r = burning rate of the solid propellant, including the 
erosive burning contribution, cm/sec 

r A = non-erosive burning rate of the solid propellant, 

0 ap , cm/sec 

Tp = flame propagation (spreading) rate, cm/sec 

T * temperature; without subscript - static gas tempera- 
ture, °K 

T a £ * average film gas temperature, (T + T p S )/ 2 » °K 

T f = adiabatic flame temperature of the solid propellant, 

1 °K 

= effective mean temperature of the igniter gas, °K 
Tpi “ initial propellant temperature, °K 
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XI 


ps 


= propellant surface temperature, °K 


{ ' 


1 

* 


s 

T 

ws 

St 

s 


t 

t. 


? . = propellant surface ^mperature, at which propellant 

PS' 1 ? ignition occurs, °K 

= duct surface temperature (in general) , °K 

= non-propellant wall surface temperature, °K 

= Stanton number 

= exponent of the length-to-diameter ratio in correla- 
tions for the heat-transfer coefficient (see Section 
D, Chapter III) 

= time, sec 

= ignition delay time (time from onset of igniter to 
the instant of first ignition of propellant) , sec 

= pressure delay time (time from onset of igniter to 
the instant of attaining 10% of the maximum pressure 
at the motor fore-end) , sec 

= time from onset of igniter to the instant of flame 
spreading completion, sec 

= time from onset of igniter to attainment of maximum 
fore-end pressure, sec 

= time from onset of igniter to attainment of m** num 
pressurization rate at motor fore-end, sec 

= gas velocity, cm/sec 

= volume, cm ^ 

= molecular weight of the combustion gases, g/g-mole 

= axial distance from igniter nozzle end, cm 

= position at the motor aft end (entrance to motor 
nozzle) , cm 

= position at the impingement of the igniter jet 

= position at the entrance to the propellant section 
(propellant slab leading edge) , cm 

pxial distance from the point of igniter jet impinge- 
ment, x - x. , cm 
un 






'fs 


'pm 


'prm 


u 

V 

W 

x 

X, 


un 


V x o 


a . ii ^ «_• 


\ 4 - 


xii 


= perpendicular distance from the propellant surface 
into the solid (see Fig. 7), cm 


Greek letters 


a 

e 

r 

Y 

A 

At 

Ax 

6 

6 

e 

e 


ph 


0 

X 

Ji 

P 

T c 

T ph 

T P 

T p,m 


thermal diffusivity, cm /sec 

erosive burning exponent [see Eq. (III-35)] 

a function of y , { y [2/ (y+1) ] / 

ratio of specific heats 

difference 

mesh size on t-axis 

mesh size on x-axis 

penetration distance (see Section G, Chapter III) , cm 
characteristic thickness of preheated layer, cm 
a small number in Eq. (III-” 1 .), °K 
equivalent sand roughness, cm 

internal performance of a solid-propellant motor 
[see Eq. (1-1)], kgf-sec/cnr 

weighting parameter for the implicit numerical scheme 

thermal conductivity, cal/cm-sec-°K 

viscosity of the combustion gases, g/cm-sec (poise) 

density (without subscript - gas density) , g/cm 3 

relaxation time for the condensed phase, sec 

time for burning through preheated layer, sec 

characteristic time for pressure variation, sec 

characteristic pressurization time, corresponding 
to maximum pressurization rate [see Eq. (C-3)], sec 


r p,t 


= total time of pressurization, sec 


characteristic time for surface reactions, sec 


xiii 


t = shear stress on the port wall, gf/cm 

w 2 
t xx = normal shear stress, g /cm 

Superscripts 

j = index for mesh time increments, t-* = jAt 


= reference quantity for non-dimensionalization (see 
Appendix A) 

= weight-averaged quantity 
= non-dimensionalized quantity 


Subscripts 

A2,A2 = values at points A ? and A_, respectively [see Figs. 
8 and 9 and Eq. (Iv-15)] 

af = evaluated at the average film gas temperature (T af ) 
d = based on diameter 

E - motor chamber aft end (entrance to motor nozzle) 
es = entrance section 
i = initial value 
max = maximum value 

n = index for mesh points in the axial direction, 

x = nAx 
n 

pd = calculated by the latest predictor calt lation [see 
Eq. IV-8) ] 

pr = solid propellant (condensed phase) 
xO = values at position x^ 

I = along right-running characteristics 

II = along left-running characteristics 

III = along particle-path line 
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CHAPTER I 
INTRODUCTION 


A. Motivation and Purpose of the Study 

Prediction and control of the pressure-time and thrust- 
time programs of high-performance solid propellant rocket 
motors during the early phases of operation often depend on 
proper understanding of the complex interactions between 
elements such as, the igniter gas flow, heat transfer to the 
propellant, flame spreading, developing flow field, and ero- 
sive burning. The ability to predict and control the motor 
starting transient enables the following typical design and 
analysis objectives to be considered and approached: 

(1) Keeping the overall transient time within speci- 
field limits; 

(2) Reducing motor-to-motor variations (which is 
important, for instance, in the start of clustered 
motors) ; 

(3) Extrapolation of sea-level motor ignition data to 
high altitude ignition conditions; 

(4) Extrapolation of tests at normal temperatures 

(- 20°c) to a wide temperature range (e.g., -40°C 
to +60°C) ; 

(5) Prediction and reduction of the overpressure, that 
sometimes accompanies ignition; 

(6) Avoiding misfires and hangfires; 

(7) Predicting the effects of long term storage on 
ignition requirements, once the effect of the in- 
fluence of aging on ignition energy is known, and 

(8) Predicting how a design modification will alter 
performance (e.g., propellant substitution, changes 
in throat area and motor dimensions, and propellant 
surface treatment) . 

Quantitative prediction and knowledge of the maximum 
pressure and pressurization rate during the ignition transient 
are very important to the rocket-motor designer. They allow 
and justify the use of a small margin of safety for the 
engine parts, thus obtaining high motor mass ratios. The 
pressure peak, frequently encountered during the starting 
transient of high-performance rocket engines, affects the 
structural design of the hardware, the grain integrity, and 
the use of acceleration-sensitive parts in the propelled 
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system. For particular designs, the rate of pressure rise 
may detrimentally affect the steadiness and ; ‘-.ability of 
burning, the viscoelastic response of the grain and inhibi- 
tors, and the dynamic response of the hardware parts. 

The broad task of ignition research is to provide a 
sound basis for defining and achieving the design and devel- 
opment goals with a minimum of expensive testing. This study 
is directed at developing analytical and experimental tech- 
niques that can be applied to achieve the above-mentioned 
goals and objectives. 

Many studies have dealt with a particular fundamental 
process taking place, during the start of solid propellant 
motors or have focused on practical correlations to facili- 
tate igniter design. Also, previous models which were 
directed at the analysis of the overall starting transient 
can be categorized into two major groups: (1) lumped chamber- 

parameter, or p(t) models, and (2) quasi-steady one- 
dimensional flow, or p(x) models. '“J- 2 The models of the 
first group assume uniform pressure and temperature distri- 
butions in the combustion chamber port. In the models of 
the second type, quasi-steady pressure, temperature, and 
velocity distributions along the port are supposed at each 
instant of time during the transient. The flame spreading 
is treated in studies of both groups in various ways: an 

experimentally determined function of time;^'7-9 partially 
or fully instantaneous; 2 #13 a linear function of the burning 
rate; ^ proceeding in a constant average rate;^ calculated 
from an experimental pressure-time plot; 3 or calculated as 
a function of the axially varying heat flux to the propellant 
surface due to changing gas temperature^ or boundary layer 
development . 5»1 5 The flame spreading analyses of Most and 
SummerfieldS and Brown et al^ 5 * are the most advanced studies 
among the lumped chamber-parameter models. A review of pre- 
vious studies is presented in Chapter II. 

Many modern, high-performance solid propellant rocket 
motors have high volumetric loading densities and small port- 
to-throat area ratios, frequently combined with large length- 
to-diameter ratios. Such motors, referred to in this study 
as HVT (High Velocity Transient) motors, are characterized by 
high internal gas velocities, significant axial pressure and 
temperature gradients, pressure overshoots, and relatively 
short transient times. The starting transients of the HVT 
motors are inaccurately predicted and unsatisfactorily ana- 
lyzed by either one of the above-mentioned classes of approach 
es. Figure 1 shows schematically a typical starting transient 
of HVT motors compared with that of motors having a high port- 
to-throat area ratio. 

This study investigates the starting transient of HVT 
motors, both analytically and experimentally. Pressure, 
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temperature, and velocity variations with time and axial 
position in the motor [p(t,x) model] are considered. Inter- 
actions between processes, such as the developing flow field, 
igniter gas flow, convective heat transfer to the solid pro- 
pellant, flame spreading, and erosive burning are taken into 
account. A practical igniter-motor combination, utilizing a 
head-end pyrogen igniter is analyzed and discussed. 

The primary objectives of the research described in this 
work, directed at achieving the previously stated tasks, were: 

(1) To develop an analytical model for calculation and 
prediction of the starting transient of HVT motors; 

(2) To design and build a laboratory-size, experimental 
HVT motor for diagnostic studies at low port-to- 
throat area ratios; 

(3) To solve the analytical model for the physical 
situation of the experimental motor; 

(4) To provide necessary empirical data and to verify 
experimentally the theoretical calculations, and 

(5) To determine the effect of different igniter, pro- 
pellant, and motor parameters on the starting 
transient. 

As can be seen from the above list of objectives, in 
addition to the many practical aspects of this work, atten- 
tion is directed at the fundamental processes that are of 
broad scientific interest, e.g., unsteady one-dimensional 
compressible fluid flows and combustion wave development. 

B. The Internal Performance of Solid Propellant Motors With 
” an Internal-Burning Grain 

Other consequences of this work are contributions to 
optimization trends in solid propellant rocketry, which are 
becoming more important as more reliable, as well as inexpen- 
sive, motors are needed for space applications (such as shuttle 
boosters and impulsive thrusters) and for sophisticated mis- 
siles. 

The internal performance of a solid propellant motor with 
an internal-burning grain, rij., may be defined as follows: 


^ Total Impulse Delivered by Motor 
i Cross-Sectional Area of Motor 



(1-1) 
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The total impulse, I 


is given by 


tot' 
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where: V is the total available volume of propellant in 

the motor , and 

I gp is the specific impulse of the popellant. 

A Volumetric Loading -Density factor, f , , is defined 
as 



(1-3) 

Substituting Eqs. (1-2) and (1-3) into Eq. (1-1) yields 


n . — p i f . l 
‘ i H pr sp vr 


(1-4) 


The first two factors in Eq. (1-4) are the so-called 
density specific impulse, Ppr^p • which describes the deliv- 
ered impulse per unit propellant volume, and is a measure of 
the quality of propellant and chemical-to-kinetic energy con- 
version process. 

The last two factors in Eq. (1-4) , f L, describe the 
amount of propellant loaded per unit cros¥-sectional area of 
the motor and are a measure of the quantity of propellant 
used, for a given motor diameter. 


This study is concerned mainly with highly-loaded motors, 
i.e., motors having high f .L-factors. High internal per- 
formance and efficient motor design require knowledge of the 
effect of the f .L- factor on the various ballistic parameters. 
The relation between propellant-charge configuration and vari- 
ous aspects of the ballistic parameters has been the subject 
of many investigations. 16 ” 19 


A concise description of this study was presented else 
where. 20 



CHAPTER II 


REVIEW OF PREVIOUS STUDIES ON THE STARTING TRANSIENT 


A large number of fundamental and applied researches have 
been conducted over the past twenty years to define and study 
the processes taking place during the ignition of solid pro- 
pellant motors and to develop methods for prediction and con- 
trol of the starting transient. Previous empirical tailoring, 
experience accumulation and rules-of-thumb are being replaced 
by analytical modeling, systematic tests, and prediction pro- 
grams. A comprehensive and orderly literature review was 
presented by Most and Summer f leld,-*'* 4 5 who summarized many 
previous studies according to approach, aspects of theoreti- 
cal model, and experimental features. Detailed data on the 
performance characteristics of typical igniter systems along 
with design criteria, as well as data on propellant response 
and prediction techniques were compiled and published by 
Falkner and Miller. 21 in the present investigation an attempt 
is made to categorize generally the studies on the starting 
transient, especially those directed at the treatment of the 
entire transient, and review briefly representative works. 

The studies on the ignition transient of solid propellant 
motors may be classified into the following four major cate- 
gories : 

(1) Studies of fundamental processes, such as ignition 
mechanisms of various propellants, chemical reac- 
tions taking place during the ignition, heat trans- 
fer to the condensed phase, flame propagation over 
the exposed surface, gas dynamics of the chamber 
pressurization, unsteady (dynamic) burning, and 
erosive burning of the propellant. 

(2) Studies of ignition systems, aimed at their optimi- 
zation by investigating the combustion properties 
of igniter ingredients, gas composition, and per- 
formance characteristics, and by formulation of 
desired compositions. 

(3) Studies of the transient of a particular rocket 
motor, when a given grain design, specified motor 
size, and/or special firing conditions (e.g., 
extreme temperatures, vacuum, or body force fields) 
require precise transient prediction and control. 

(4) Studies of the entire ignition transient, directed 
at the calculation of the pressure-time and thrust- 
time histories, considering interactions between the 

previously described transient processes and simul- 
taneous phenomena. 
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Some studies belong to two or more categories. Another 
classification of studies may be made according to the labora- 
tory-type of energy stimulus, which determines the mode of 
energy transfer to the propellant: shock tube, arc image 

furnace, radiation furnace with hot stagnant gas, laser appa- 
ratus, pyrotechnic or pyrocren-type igniter, operating in a 
real rocket motor, etc. This study belongs to Category (4). 

A. Studies of Fundamental ProceF-ses 

Most previous studies of the starting transient belong 
to the first, above-mentioned category. Comprehensive criti- 
cal reviews of contemporary ignition theories are given in 
Refs. 22, 23, and 24. Analytical and experimental studies 
of the mechanisms of solid propellant ignition have resulted 
in three major ignition models, that are classified as solid 
phase, gas phase, and heterogeneous ignition theories. These 
theories differ mainly in the postulated place of the domin- 
ant, exothermic, runaway reactions, physical state of reac- 
tants, and sequence of processes. Accordingly, different 
ignition criteria may be presumed. The solid phase ignition 
theory, analytically established by Frazer and Hicks. 23 ' 26 
and further developed at the University of Utah, 2 ^ -2 ® an( j a t 
NWC (formerly NOTS) 3 ®' 33 assumes, that condensed phase exo- 
thermic reactions dominate the runaway process, and neglects 
any participation of gas phase species. The gas phase igni- 
tion theory, developed mainly at Princeton University , 32-34 
assumes that fuel vapors created from decomposition-gasifica- 
tion of the solid fuel diffuse into the hot, oxidizing, 
environmental gas and react exothermal ly near the propellant 
surface. -The heterogeneous ignition theory, developed mainly 
at UTC 33 ' 3 ® and recently advanced by Bradley and Williams, 7 ' 3 ® 
and by Waldman and Summerf ieid 3 ® t 40 proposes, that the hot 
oxidizer gas formed from the decomposition of the ammonium 
perchlorate oxidizer reacts exothermally with the fuel in the 
condensed phase to provide the controlling reaction. The 
chemical reactions, involved in the ignition and combustion 
processes, have been treated in a global manner, 2 ®' 3 ®' 33 or 
in a general scheme, according to the model, 43 or by extensive 
investigation of a certain process, such as the decomposition 
of the ammonium perchlorate oxidizer, 42 *" 44 and other propel- 
lant ingredients. Experimental studies try to determine gas 
species composition and concentrations by spectrometry and 
mass spectroscopy. 

The process of heat transfer to the solid propellant, 
in all its major modes, has received considerable attention, 
due to its vital importance to the ignition transient. The 
rate of propellant heat-up is a complex function of simultan- 
eous heat-transfer processes, such as convection, radiation, 
hypergolic heating condensation, hot-particle impingement, 
and chemical recombination. Moreover, the unsteady turbulent 
flow with chemical reaction over the rough, decomposing pro- 
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pellant surface, in a complicated-geometry duct, makes an 
appropriate analytical study a formidable task. Instead, 
current analyses use simplified, semi-empirical heat-transfer 
correlations, based on experimental results. In the follow- 
ing paragraphs, numerous published studies are described 
briefly. Some heat transfer phenomena and correlations are 
discussed further in Chapter III. 

An extensive research program, studying various aspects 
of ignition in practical situations has been conducted at the 
United Technology Center since 1962. Pyrogen igniters, both 
head-end and aft-end, have been fired into instrumented solid 
propellant motors, 45 and into a cylindrical copper duct, 
equipped with an exhaust orifice. ^6 Heat-flux and flame- 
propagation measurements were made in a series of experiments, 
which resulted in empirical heat-transfer correlations for 
both the induction and flame spreading periods. ,A1 These 
correlations are based on expressions for fully-developed 
turbulent flow, modified for duct-entrance and jet-impinge- 
ment effects. Analytical models for ignition and flame 
propagation, taking into account contributions of exothermic 
heterogeneous chemical reactions and assuming propellant 
surface autoignition temperature also were established.^ ,4 5 
The theory considers axial variation of the heat flux due to 
the drop in gas temperature along the port and increasing 
convective fluxes during flame spreading. Agreement between 
theoretically predicted and experimentally measured ignition 
delay times and flame propagation rates were reported. Ex- 
tensive data on the effect of heat-flux, pressure levels, 
and solid-phase composition on the ignition response of vari- 
ous propellants also have been published by UTC.^®'^® 

Carlson, Seader, and Wrubel of Rocketdyne carried out 
a thorough study of the heat transfer from pyrogen igniters 
to propellant grains of various configurations. First, they 
used a simulated igniter-motor combination, in which a hot 
nitrogen gas was forced through an instrumented duct. 50-52 
The duct consisted of thin copper liners (installed in stain- 
less steel test sections) , and had a nozzle at the aft end. 
Various head-end and aft-end igniter-nozzle configurations 
were tested. Later, practical pyrogen igniters were fired 
into the instrumented duct and into real propellant charges. 53 
Measured heat-transfer coefficients were compared with the 
classical turbulent pipe flow correlation and plotted versus 
nondimensional distance. The results are of particular 
interest to this study, since port-to-throat area ratios of 
1.1, 2.0, and 3.0 were used. Unfortunately, no analytical 
correlations have been deduced by the investigators. However, 
their research revealed many of the parameters affecting the 
heat transfer to propellant grains, such as, entrance effects, 
igniter jet structure, port configuration and port-to-throat 
area ratio. 
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Lovine and Fong, at Aerojet-General Corp. have also 
measured heat fluxes by firing subscale, pyrotechnic igniters 
into an instrumented test chamber. They correlated their 
heat-transfer data for the inlet region of the chamber by the 
expression for turbulent flov over a flat plate. However, a 
space-averaged heat-loss term is used in their gas-dynamics 
equations. Later, Micheli and Linfor^S developed separate 
expressions for the heat-transfer coefficient for the three 
distinct regions of the port i ' a real rocket motor, namely, 
the impingement, downstream-of-impingement , and u t ~-ream-of- 
impingement zones, including a radiation contribution. 

The ignition of solid propellants by radiative, convec- 
tive, and chemical heating was the subject of many studies 
at Arthur D. Little, Inc. 56-59 Convective heating was stud- 
ied by two methods: (1) in a hot-gas tunnel, by passing com- 

bustion gases through a rectangular test section, in which a 
propellant sample was held; and (2) in a cylindrical instru- 
mented duct, into which practical pyrotechnic and pyrogen 
igniters were fired. Heat-flux measurements were interpreted 
in terms of developing turbulent flow correlations, assuming 
steady-state flow of the igniter ras. A simple two-layer 
film model, which accounts for the combined effect of all 
non- radiative heat-transfer modes, was '.ntroduced. Effects 
of gas velocity and port configuration were qualitatively 
studied by Bastress and Niessen.59 They observed a marked 
increase in ignition times the flow velocity is increased 
from Mach number 0.29 to 0.46. For long propellant samples 
first ignition was observed downstream from the sample lead- 
ing edge. From the accumulated data Bastress and Niessen 
concluded that ignition behavior is sensitive to the local 
concentration of species evolved from upstream portions of 
the sample. 

Intensive studies on radiative and convective ignition 
of solid propellants have been carried out at Princeton Uni- 
versity. McAlevy 32 and Kashiwagi et al 6 ^ studied the igni- 
tion of specimens of a solid fuel or propellant exposed to 
a hot, high-pressure stagnant oxidizing gas in a shock-tube. 
Their findings support the gas-phase ignition theory, and can 
be used as input to an applied ignition theory suitable for 
analyzing practical motors ignited by hot non-oxidizing high- 
velocity gases. The heat flux in an experimental rocket 
motor using a hot-wire, or gaseous pyrogen igniter was measur- 
ed by Parker and Most . 3 • 15 • 61 • ® 2 A correlation of the type 
Nux * a^Re x 0*R was deduced from two measurement positions, 
large scatter of data, and several ad-hoc assumptions, such 
as, average constant driving temperature difference and a 
leading edge parameter based on a dif f icult-to-def ine duct 
leading edge. Both investigators report, however, good agree' 
ment between measured and predicted (with the above-mentioned 
correlation) flame spreading rntes. 



Extensive and systematic studies of AP-propellant igni- 
tion processes, by radiative heating in a radiation furnace, 
and convective heating in a shock-tube test section, have 
been conducted at the University of Utah. Baer, et a l 27 ' 28 
measured convective heat fluxes and ignition times of propel- 
lant samples, mounted in the wall of a test section at the 
end of a shock tube. They found that ignition time depends 
primarily on propellant composition and heat flux, when air 
or nitrogen was used as the test gas. In this case, the 
dependence of ignition time on heat flux was correlated ac- 
cording to the solid-phase thermal ignition theory. When 
the propellant was heated by oxygen, a significant reduction 
in ignition delay, attributed to exothermic reaction between 
the 9J$y?g n and the fuel, was observed. Keller, Baer, and 
Ryan® 3,6 in a systematic experimental study, investigated 
the effect of gas velocity on ignition, discovered previously 
by McCune 2 ® and Kling et al.^5 High-pressure (14 to 25 atm), 
high-temperature (I000-2600°K) gas behind the reflected shock 
wave in a shock-tube end was passed through a constant-area 
rectangular test section. The propellant samples were mount- 
ed flush with one of the test-section walls, and the test-gas 
velocity was controlled by critical-flow orifices at the down- 
stream end of the section. Test-section Mach number was 
varied over the range of 0.07 to approximately 1.0. In gen- 
eral, the effect of higher gas velocities was to decrease 
the time to ignition for the same heat flux. This effect 
was found to decrease at very high velocities. Keller, Baer, 
and Ryan discovered an interesting effect of the propellant- 
surface condition, i.e., for the same heat flux at the sur- 
face, reduced ignition times were observed, when propellant 
surface was rougher, and gas velocity was lower. For pro- 
pellants with a smooth surface, the effect of gas velocity 
on ignition was much diminished. For propellants with a 
rough surface, the effect of gas velocity was substantial 
only for gas velocities under Mach 0.5. It was postulated, 
that the improved ignitability of propellants with rough 
surfaces is produced by supplementary secondary ignition 
reactions among propellant decomposition products, that are 
generated at the surface by increased heating of surface ir- 
regularities. If the velocity of the heating gas is increased 
the decomposition products are swept into the heating gas 
before reacting. It is unfortunate, that pressure effects 
were not considered in their study and no experiments were 
carried out at high velocities and heat fluxes, and low pres- 
sures, pertinent to HVT motors. In later studies by Richard- 
son, Ryan, and Baer, 66, 67 the surface roughness effect was 
not observed at lower heat fluxes (2 to 30 cal/cm 2 -sec) and 
lower test-gas temperatures and pressures (760 to 1300°C, 2.9 
to 7.7 atm, Mach Number from 0.05 to 0.3). Also, the obser- 
vation of Bastress and Niessen^S that the first evidence of 
ignition occurred at some distance from the leading edge was 
not confirmed, under these test conditions. The effect of 
test-gas temperature was also studied. Convective heat fluxes 
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were measured in many studies at the University of Utah ' ° ' 
66-68 a nd the heat-transfer coefficient was correlated in the 
form h c = const * T0*3G m . The power m of the mass velocity 
G was reported to have values ranging from 0.68 to 1.06, in- 
dicating, that in some experiments transition flow may have 
existed over the propellant sample. The values of the con- 
stant and the power m were found to depend upon the nature 
of hot gas used in the tests. It should be noted, that un- 
steadiness of the convective heat transfer during tests under 
constant flow conditions (7.7 atm pressure, Mach Number of 
0.08, gas temperatures of 760 and 1000°C) was reported by 
Baer and Ryan. 67 This unsteadiness, which was not observed 
with different measuring techniques and at higher gas temper- 
atures, v.*as attributed to nonuniform flow in the test section. 
In another study of the effect of propellant surface rough- 
ness^ propellant samples with surface imperfections as large 
as 30 microns were subjected to convective heating in a shock 
tube and the heating was stopped before ignition. These 
same samples were then immediately subjected to ignition 
tests. The results did not show agreement with the hypothe- 
sis that surface roughness causes early ignition, but no 
definite conclusion has been drawn. 

An investigation of the efffct of environmental gases 
and pressure on the ignition of u solid rocket propellant 
and ingredients was conducted by Hightower. 70 The propellant 
samples were ignited in an arc-image furnace by radiant energy 
fluxes over the range of 6.5 to 120 cal/cm^-sec. Environmen- 
tal gases used were air, oxygen, nitrogen, methane, and 
helium. The tests were carried out at pressures of 15, 45 
and 100 psia. At the higher energy fluxes the ignition time 
was found to be dependent on the environmental pressure with 
ignition times increasing when decreasing pressure. The 
shortest ignition times were obtained with pure oxygen as the 
environmental gas. Hightower concluded that, under the ex- 
perimental conditions of his investigation, exothermic gas 
phase reactions occurring in a zone immediately adjacent to 
the propellant surface play a role of major importance in the 
ignition of AP composite propellants by high intensity ther- 
mal radiation. 

The rate of flame spreading over the exposed propellant 
surface is a necessary element for any analysis of the start- 
ing transient. Flame spreading is a complicated process, 
coupling igniter behavior, port geometry, propellant response, 
unsteady combustion, and chamber gas dynamics. It has been 
studied both experimentally and analytically, coupled, or un- 
coupled with other processes. 

A comprehensive theoretical model of flame propagation 
along a smooth port was developed at UTC. 1 '*'^ This model 
takes into account all major possible modes of heat transfer 
to the propellant (i.e., conduction, convection, radiation, 
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and heterogeneous exothermic reactions) , and considers tem- 
poral and axial variations of the gas temperature and heat 
transfer. During the flame spreading the convective heat- 
transfer coefficient is enhanced by an experimentally-deter- 
mined function of the igniter and burning propellant mass 
flow rates, h c = h C/ig [(m^ g + p pr rA b ) /mi g ] 0 , where h ? ig 
is the heat-transfer coefficient for igniter operation al6nfc, 
correlated for developing^ or developed^ turbulent flow. 

The flame spreading process is determined by calculating the 
time to reach a critical surface temperature for ignition 
for positions along the longitudinal axis. Besides use of 
empirical factors to account for the enhancement of the heat- 
transfer coefficient, this model does not consider coupling 
to the transient gas dynamics of the chamber. Instead, 
neither temporal nor spatial variation of the chamber temper- 
ature were considered, and only temporal change of pressure 
was assumed. Fair to good agreement between calculated and 
measured (by the trip-wire technigue) flame spreading rates 
was reported. ^ 5 in another study at UTC,"^ under the un- 
realistic assumptions of no axial variation of the heat-trans- 
fer coefficient and propellant surface temperature, a scaling 

1 A ^ 4— W a P a -w «• A . /*» ^ P 4- / t.rVi a A v* 4 a* 4- V* 4 


law of the form 
spreading rate. 


q£d p /St 


(where 


r p is the flame 


is Stanton Number, and 
port) was derived. 


qg, is the heat flux to the propellant, St 


is the diameter of a cylindrical 


2 72 

DeSoto and Friedman ' developed in 1964 a model for 
flame spreading, which assumes an isothermal one-dimensional 
gas flow in the chamber, spatially uniform pressure and mass 
generation due to the grain burning only. The initial con- 
ditions in their analysis include instantaneous ignition of 
part of the grain, uniform pressure and temperature in the 
gas phase, a uniformly distributed temperature in the unig- 
nited portion of the propellant surface, and an axially 
varying velocity distribution. Pressure varies with time 
only; velocity varies with time and axial distance. The 
analysis was used also for calculation of the chamber-filling 
period. No comparisons with experiments were reported. In 
view of the many unrealistic assumptions, the practical appli- 
cation of this analysis is doubtful. 


13 

Paul, Lovine, and Fong, of Aerojet-General Corp. , pre- 
sented an analysis in 1964, in which the area subjected to 
direct igniter jet impingement is assumed to be ignited 
instantaneously. The heat transfer downstream of the impinge- 
ment zone is taken as being independent of axial distance, 
the gas temperature is assumed to be constant with respect 
to time and space, and only the temporal variations of the 
uniform chamber pressure and heat flux are considered. Adopt- 
ing the critical-temperature ignition criterion, this analy- 
sis predicts the time to instantaneous flame spreading down- 
stream. However, in practice, instantaneous flame spreading 
is generally not observed. The flame spreading in the stagnant 
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region upstream of the igniter jet impingement zone is cal- 
culated by a simplified analysis, assuming temporal and 
axial variation of the convective heat transfer and consider- 
ing . adiation from a cloud and heat conduction, as well. 

Considerable research on the flame spreading and chamner 
filling processes has been conducted at Princeton University. 
A model, established by Parker, Summerf ield , et al, 3, 61, 73-7 5 
is based on the so-called dynamic temperature, well-stirred 
analysis. This analysis assumes uniform chamber pressure 
and temperature, both changing with time. The first method 
of using the analysis uses a given experimental pressure 
variation and an assumed initial temperature. A system of 
two ordinary differential equations for the time variation 
of pressure and temperature is solved numerically to yield 
the flame- spreading rate. In a second method of using the 
analysis, initial conditions of burning area, chamber tem- 
perature and propellant surface temperature are assumed. 

In addition, an experimentally determined heat-transfer cor- 
relation and a critical temperature for ignition are used. 

The method solves the pressure and temperature variation 
equations and the heat equation for the solid phase to yield 
the chamber pressure history and flame spreading rate. Con- 
ditions of gasless igniter are assumed in Parker's analysis. 
The aforementioned assumptions of initial conditions make 
the use of this analysis for prediction of the flame spread- 
ing and pressure transient rather inconvenient. Parker's 
analysis was modified by Most, Summerfield, et al . #15,76-78 
Most's analysis includes the igniter operation until an early 
stage of the flame spreading (rather unrealistic assumption) , 
and calculates the flame spreading rate by solving the solid 
phase heat-up equation, coupled to the gas-dynamics equations. 
The heat flux is both time- und distance-dependent. However, 
the driving temperature difference is assumed to be constant 
in both time and distance. The ignition criterion used in 
Most's study was also that of a critical surface ignition 
temperature . 

The flame spreading over the surface of double base 
propellants and propellant ingredients in a quiescent atmos- 
phere was studied by McAlevy, et al.^°'®^ It was found ex- 
perimentally that flame-spreading velocity varies directly 
with pressure (between 0.1 and 1 atm), oxygen fraction of 
environmental atmosphere, and surface roughness. The flame 
spreading was presented an a continuous gas-phase ignition 
process and the analytical model derived was verified exper- 
imentally for smooth surfaces. For rough- surf aced samples, 
photographic evidence of random ignition sites, ahead of the 
flame front, was obtained. 

81-8 3 

Mitchell and Ryan conducted flame spreading experi- 

ments with a rarefaction-tube apparatus, under conditions of 
constant pressure and velocity using a relatively cool gas 
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flowing over the propellant-sample surface. The T-shaped 
propellant sample was positioned near the nozzle-end of the 
tube and flame spreading was observed by means of high-speed 
photography. In their analysis, the effects of pressure and 
velocity, both enhancing the flame spreading rate, are isola- 
ted. The effect of pressure was found to be greater than the 
effect of velocity for a given fractional change. 

In other studies, as mentioned in the Introduction , the 
flame spreading rate is based on very unrealistic assumptions, 
such as having an average constant value, ® or being a linear 
function of the burning rate. In some investigations, flame 
spreading rates are scaled from experimental data for similar 
systems . 1 

Recent studies on the response of solid propel la'~ 4_ burn- 
ing to pressure changes have shown that during rapid pressure 
excursions (pressurization or depressurization) the instantan- 
eous burning rate differs from the ste ady-stat ' values, cor- 
responding to the instantaneous values of pressure. During 
the pressurization process in the starting transient, dynamic 
burning rates higher than the steady-state may be encountered 
due to unsteady thermal profile in the solid phase. Early 
investigations of the effect of pressure transients on the 
burning rate resulted in Von Elbe's model and formula for the 
dynamic burning rate. 8 ^/ 85 This model and formula were sub- 
sequently modified and used by Paul, Lovine and Fong®® at 
Aerojet and others. Detailed models of the combustion wave 
for calculation of the dynamic burning rate response were 
developed by Merkle, Turk, and Summer field, 8 ^ by Krier, T’ien, 
Sirignano, and Summerfield (KTSS model) , 8 ® by Denison und 
Daun, 88 by McClure, Hart, and Bird, 88 etc. A critical class- 
ification and review of the various analyses was presented 
by Culick. 81 Recently, an approach, established by Zeldovi.ch 
and extended by Novozhilov, 92-94 W as successfully applied at 
Princeton University by Summerfield, Caveny, Battista, et al 95 
97 for calculating extinction boundaries, and transient inter- 
nal ballistics. In this latter method, experimental ste^ady- 
state burning rate data [r(p,T ^)] over a wide range of 
initial (ambient) temperatures^is needed 88 ' 8 ^ instead of a 
detailed flame structure. It has been shown 8 ® • 88 • 8 ? that 
the dynamic burning causes pressure overshoots during igni- 
tion. 


Another aspect of dynamic burning is the burning of pre- 
heated propellant layers during the flame spreading and cham- 
ber filling periods.®! This is especially pronounced 
after long induction (heat-up-to-ignition) periods and may 
cause greatly increased instantaneous burning rates and sub- 
sequent pressure overshoots. 

Erosive burning is a phenomenon frequently encountered 
in HVT motors. This process of burning rate increase, due to 


the flow of combustion gases over the propellant surface, is 
most pronounced during the chamber filling period. Like the 
flame-spreading rate, quantitative knowledge of the erosive 
burning is necessary for the prediction and control of the 
starting transient. Therefore, numerous studies have been 
conducted on the subject, and many empirical and Femiempiri- 
cal correlations were developed in attempts to accurately 
account for erosive burning. Unfortunately, a totally- 
successful and general model of erosive burning has yet to 
be developed. An extensive research program on erosive burn- 
ing was carried out at Purdue University by Zucrow, Osborn 
and Murphy . 98-101 ^ comprehensive review of previous inves- 

tigations and laboratory methods was prepared by Murphy, 1 Gl- 
and by Williams, Barrere, and Huang . 1 '-* 2 The general conclu- 
sions drawn from the experimental studies are: ( 1 ) slow- 

burning propellants are more prone.to erosive burning than 
fast-burning propellants ;*®' 1 ' ( 2 ) there is a threshold 

velocity (or mass velocity) for erosion, which may be pres- 
sure-dependent and is more pronounced for the faster burning 
propellants ; 1 04 -1( n ( 3 ) arain configuration affects the 

erosive burning; 10 ^, 107,108 ( 4 ) at relatively high pressures 

(100 atm) the type of propellant binder and the binder-to- 
oxidizer ratio have little effect on the erosive burning, 
but at medium and low pressures (50 and 20 atm) the effect of 
binder type is large j 1 ^, 109 and (5) addition of aluminum 
has only a slight effect on the erosive burning sensitivity 
of o pro f ellants based on polyurethane and UTREZ binders. 108, 

Among the analytical approaches, ' 112-11 5 that of 
Lenoir and Robillard 112 is the most widely accepted and 
experimentally confirmed. 12 ,105-107 ,111,116 The theory of 
Lenoir and Robillard postulates that the erosive-burning 
increment of the burning rate is proportional to the coeffi- 
cient of convective heat transfer from the core of flowing 
combustion gas to the propellant burning surface. The latter 
is attenuated by the blowing from the burning surface accord- 
ing to an exponential law, derived from the data of Mickley, 
et al 11 ' and Rubesin, et al, 11 ® and thus the total burning 
rate law takes the form r = ap n + kh c exp(- 8 rpp r /up) , where 
k is the proportionality coefficient evaluated analytically 
and 8 is an experimentally determined constant. Lenoir and 
Robillard used the Chilton-Colburn correlation for the heat- 
transfer coefficient, valid for turbulent flow over a flat 
plate. Lawrence, et al 111 used successfully Lenoir-Robillard ' 
law with a turbulent pipe-flow heat-transfer correlation for 
grain length to port diameter ratios of up to 50. 
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B. Studies of Ignition Systems 

Many studies on the igniter hardware and characteriza- 
tion of igniter materials for the initiation and energy re- 
lease systems have been conducted in the induscrial and 
governmental laboratories. Ignition systems fall into one 
of the three general types: pyrotechnic, pyrogen, and 

hypergolic. Each type has its own features and way of devel- 
opment. All three categories of systems, as well as multiple- 
ignition systems for use with restartable or pulse solid- 
propellant motors, are the subject of past and present exten- 
sive activity. A comprehensive review of the work, done by 
the major agencies, engaged in research and development of 
ignition systems prior to 1956, was reported in a concise 
form by Korpi. 1° Up to 1956, a great deal of practical work 
on development and performance of pyrotechnic systems, and 
formulation of igniter materials, such as metal-oxidizer 
compositions, was done in Aerojet, 120 Thiokol, Alleghany 
Ballistics Laboratory (ABL) , Picatinny Arsenal, Ballistic 
Research Laboratory (BRL) , Jet Propulsion Laboratory (JPL) , 
Naval Ordnance Laboratory (NOL) , Naval Ordnance Test Station 
(NOTS) , University of Michigan, University of Utah, etc. In 
the past fifteen years much development work has been con- 
ducted by the Bermite Powder Co. (BPC) . 

A monograph on solid-propellant rocket igniters recently 
published by NASA, -*-21 presents accumulated experience and 
knowledge in an organized concise form for effective use in 
design. A highly recommended igniter-design handbook ! 22 
which compiles experience gained and practical information, 
was prepared by Bermite Powder Co. (BPC) . 

Many studies have tried to determine the minimum igniter 
energy, required for propellant ignition, which, in turn, 
determines the igniter system characteristics. In the Bureau 
of Mines, ignition energy values from 2 to 15 cal/cm 2 for 
double base propellants were found. 19 *123 studies at BPC 
and NOL have resulted in empirical relationships between 
motor configuration parameters, propellant ignitability , type 
of igniter material used and energy required for satisfactory 
ignition. 5 > 121 pointed out by several investigators (e.g., 
Refs. 5, 23-27, 63, 64, 66-69), the ignition limit is estab- 
lished by the rate of ignition energy supplied, rather than 
by the total energy only. Other studies on igniter sizing 
are based on attaining a given necessary-for-ignition pressure 
in the motor port. 124 

Studies of the energy transfer supplied by the igniter 
gases by convection and radiation (the latter i very pro- 
nounced for usual pyrotechnic compositions) were conducted 
by firing the igniters, both pyrotechnic and pyrogen-type of 
different charge weight, into instrumented ducts that simu- 
late solid-propellant motors . 125 , 126 
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A brief listing of igniter material considerations and 
applications, as developed at BPC, was recently presented by 
Robertson. 1^7 Hypergolic ignition systems, applying a hyper- 
golic reagent (usually liquid chlorine trifluoride) to the 
propellant surface were developed mainly at UTC.128,129 

C. Studies on the Transient of Particular Motors 

The effects of motor and igniter configuration, igniter 
placement, and initial conditions on many processes taking 
place during the ignition transient, make special studies of 
particular cases inevitable. For instance, a specific grain 

design may cause the flow pattern of the igniter and propel- 

lant combustion gases to produce transients differing signif- 
icantly from those obtained using a grain with a cylindrical 
port of the same cross-sectional area. Studyl^O j_ n which 
pyrogen igniters were fired into instrumented star-shaped 
port has shown that the heat flux at the star inner diameter 

is about 3 times higher than that at the outer diameter. 

Many transient prediction and design programs were developed 
for specific grain design and motor configuration. This is 
why every correlation model and analysis has to be checked 
for its domain of applicability, before being applied to a 
given situation (see p. 18 in Ref. 5) . 

The development of the Algol III solid-propellant 
motorUG for example, was influenced by the configuration 
dependence of the erosive burning, which affects the starting 
transmit. Inhibiting part of the exposed propellant surface, 
in this case, has brought the starting transient into an 
acceptable range. Conditions of spin and altitude have 
affected the performance and the transient of many rocket 
motors . 131-133 For th e same igniter mass flow rate, motor 
spin shortens the transient and increases considerably the 
pressure level. The existence of a nozzle plug or diaphragm 
changes the entire pattern of flow in the chamber and heat 
transfer to the propellant, affecting largely the pressure- 
time curve. ^ High longitudinal motor acceleration may affect 
considerably the transient motor internal ballistics by reduc- 
ing the pressure drop along the port as well as the fore-end 
pressure. These are just a few examples to show the effects 
of particular conditions, which make the general application 
of an even comprehensive analysis a difficult problem. 

D. Studies on the Entire Ignition Transient 


As mentioned in the Introduction, previous studies on 
the entire starting transient may be divided into two major 
classes: (1) the lumped chamber-parameter models, known also 

as the homogeneous models, and (2) the quasi-steady, one- 
dimensional flow models. The major difference between the two 
groups lies in the assumptions under which the analytical model 
of the chamber gas dynamics is formulated. 
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The analysis of Adarns^"' presented in 1966 is typical 
of the early studies of the first aforementioned group. The 
gas temperature in the chamber is assumed to be uniform space- 
wise and constant timewise (the so-called "isothermal analy- 
sis") . The pressure transient is calculated from the mass 
conservation equation for the motor chamber, assuming uniform 
pressure, and caking into account the igniter mass inflow 
rate. Flame-spreading rates are presumed to be a prescribed 
function of time (usually a linear function) and ignition 
times are also provided. Erosive burning effect is accounted 
for in the analysis by taking it as a function of an average 
value of the Mach number along the port axis. Burning rate 
correction for pressurization rate is included. Experimental 
values of c* are used to account for heat transfer to the 
propellant. Reasonable agreement between calculated and 
experimental pressure-time histories is reported and it is 
claimed that ignition peaks are caused primarily by erosive 
burning. All assumptions in this analysis are inadequate for 
HVT motors. Mainly, the assumptions of uniform pressure and 
constant chamber temperature and the concept of an average 
Mach number are unsatisfactory. 

Another isothermal analysis was presented by DeSoto and 
Friedman^/72 i n 1964 . in their analysis, however, the flame 
spreading rate is calculated by consideration of spatial, as 
well as, temporal variation of the gas velocity and heat 
transfer to the propellant, and the use of the critical sur- 
face-temperature concept as ignition criterion. Thus a con- 
tradiction between the isothermal assumption, and the varying 
heat transfer is created. No igniter mass flow rate is con- 
sidered and no comparison between calculated and experimental 
pressure transients is reported. 

136 

Sharn, et al (1964) carried out an experimental and 
analytical study of the ignition transient by firing a pyro- 
gen igniter into a small motor. The ignition time is calcu- 
lated using experimental data and the concept of critical 
surface temperature for ignition. An isothermal analysis, 
combined with an instantaneous flame spreading is used for 
the calculation of the pressure transient. A single compar- 
ison with an experimental pressure-time history shows good 
agreement during the first portion of the transient and a 
poor agreement during the later part. 

An isothermal analysis of the starting transient was 
developed by Brown, et al 1 '*,45 a t UTC. The analysis uses 
time-dependent burning surface area which is determined from 
a flame-spreading model. The latter considers axial varia- 
tion of the gas temperature and heat flux to the propellant. 

A further development of the lumped chamber -parameter 
models was the use of uniform but timewise changing chamber 
gas temperature (the so-called "dynamic temperature analysis") . 
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This is performed by combining the mass and energy conserva- 
tion equations for the gas in the combustion chamber. The 
’’dynamic temperature analysis" was developed largely at UTC 5 
starting with the analysis of Baker 137 in 1962. However, no 
attempt to couple the more-complicated chamber gas dynamics 
to a calculation of the flame spreading rate was reported. 
Instead, experimentally obtained time to first ignition and 
flame spreading rates or a specified burning surface area as 
a function of time were used. 137 The analysis of Baker con- 
siders the regimes of both choked and unchoked motor nozzle. 

It can also be used to calculate flame spreading rates for 
a given chamber pressure transient. Baker emphasized the 
error introduced by assuming constant stagnation temperature 
in the motor. He used his analysis to predict the starting 
pressure transient of a slotted-grain motor and reports good 
comparison with measured transients. 

Bradley 4 published in 1964 a dynamic temperature analy- 
sis which assumes a homogeneous, adiabatic system. This was 
the first analysis to predict the effect of various igniter 
and motor parameters, e.g., burning rate pressure exponents 
of the igniter and propellant, igniter flame temperature, 
igniter-to-propellant mass flow ratio and the initial pres- 
sure and temperature in the motor. A general igniter form 
function was used to show the effect of varying igniter mass 
flow rate on the transient. Assuming a linear relationship 
between the flame spreading rate and the burning rate, 

Bradley showed the importance of flame spreading in the 
pressurization process. In particular, he pointed out the 
large effect of the fraction of propellant surface initially 
ignited. 

54 

Lovine and Fong used a dynamic temperature analysis 
fcr the chamber filling period only, combined with expressions 
for unsteady burning rate and unsteady nozzle flow. The in- 
duction interval was treated by evaluating average values of 
igniter-induced heat flux and port pressurization and by using 
experimental pressure-heat flux-ignition time data. The 
flame- spreading interval was not treated in Ref. 54. In an 
earlier section it is mentioned that in another report 13 Paul, 
Lovine and Fong consider an instantaneous flame spreading. 

The pressure transient analysis established at Princeton 
by Parker and Summer field 3 / 61,7 3 and later modified by Most, 
Summerfield, et al^' 15 ' 7 ®' 7 ® can be considered as the more 
comprehensive of homogeneous models. Parker's analysis, con- 
sidering only the flame spreading and chamber filling inter- 
vals is described earlier in this Chapter. Since Most's work 
was the precursor of this study it deserves more description 
and attention. Most's analysis calculates the entire ignition 
transient. The three major assumptions on which his model is 
based are: (1) uniform pressure in the chamber (Vp * 0) ; 

(2) uniform gas temperature in the port (VT = 0) , and (3) 
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kinetic energy of the combustion gases is negligible, com- 
pared with their thermal energy (u 2 << c p T) . All three 
assumptions are inappropriate for HVT motors (having low 
port-to-throat area ratio) . Assumption (2) was found to be 
unrealistic by many investigators, including Most, even for 
motor with large A p /A t due to the high heat-transfer rate 
in the port, causes by large temperature difference. Most 
accounts approximately for this by using a lower constant 
igniter gas temperature, determined experimentally. In 
Most's analysis the two nonlinear ordinary differential equa- 
tions for the chamber pressure and temperature are coupled 
to the flame spreading process through the instantaneous burn- 
ing area. The latter is calculated from the solid phase heat- 
up equation, using an experimentally determined heat-transfer 
coefficient, and the critical ignition temperature, as the 
ignition criterion. The theoretical model was tested exper- 
imentally over a wide range of parameters. Good agreement 
between the calculated and measured pressure transients is 
reported. 

15 

Most reported that Falkner and Kilgroe of CETEC 
Corp.l 2 ** had used his model with reasonable success, utiliz- 
ing heat-transfer correlations of the form developed by 
Kilgroe. 47 

Motors with high volumetric loading density and low port- 
to-throat area ratio have been analyzed by the one -dimensional 
quasi-steady models. In these models, a series of con- 

trol volume increments is assumed along the port axis, and 
one-dimensional (steady-state flow) conservation equations 
are written for the increment (the so-called "incremental 
analysis"). Variations of the port area, burning surface and 
burning rate along the port are taken into account. 

7—9 139 

The usual calculation technique ' in the quasi- 
steady models starts with an approximation of the fore-end 
pressure from the steady-state formula for equixibrium pres- 
sure. Then, the steady-state continuity and momentum equa- 
tions are iterated for each increment using the perfect gas 
law and the steady-state energy equation to yield the flow 
parameters at the exit plane of the increment. Nozzle en- 
trance stagnation pressure predicted by this method is com- 
pared with that necessary for passing the calculated total 
mass flow rate through a given nozzle-throat area. The iter- 
ation procedure repeats until the two stagnation-pressure 
values converge. This calculation is repeated for each time 
step in the transient. All calculations are organized in 
computer programs, the descriptions of which are given in the 
aforementioned references. No coupled calculation of the 
flame spreading rate has been given by any of the reported 
analyses. The flame spreading is required as a prescribed 
input to the computer calculations. The induction interval 
is not treateu. The computer program in Ref. 8 can calculate 
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ignition-delay times if the spatial and temporal distribu- 
tion of the heat flux created by the igniter is specified. 

No detailed comparison of calculated and experimental pressure 
transient histories is reported. By their character, the 
quasi-steady models attribute pressure peaks during ignition 
only to erosive burning and mass acceleration losses along 
the port; for they contain no provision to account for the 
dynamic effects which may be prominent in KVT motors. 

In Ref. 10 an attempt is made to account for unsteady 
gas dynamics in the port. The analysis was aimed at the pre- 
diction of the starting transient of the HIBEX motor, charac- 
terized by high axial acceleration forces and use of fast 
burning propellant. Unsteady mass and momentum conservation 
equations with spacewise variation of pressure, temperature 
and velocity, combined with the perfect gas law, are derived 
in the incremental analysis. However, in the solution by 
finite-difference methods, the timewise derivatives of tem- 
perature and velocity are neglected and use is made of the 
quasi-steady energy conservation equation. Good agreement 
between calculated overall fore-end pressure-time plot and 
three experimentally measured pressure histories is presented, 
but there is no detailed comparison of the starting transient 
part of the plots. 
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CHAPTER III 

ANALYSIS OF THE STARTING TRANSIENT OF HVT MOTORS 


A. Description of the Physical Model 

The analytical model is established for the general 
class of HVT motors and is solved for the physical situation 
corresponding to the experimental rocket motor designed and 
built (as part of the study) to provide detailed information 
on the ignition events and to test the validity of the ana- 
lytical model . 

The physical situation, analyzed and solved by the model, 
considers a rectangular window motor, which consists of an 
entrance section, a propellant section and an aft section, as 
described in Fig. 2. Two opposing propellant slabs (2.54 cm 
wide, 0.635 or 0.318 cm thick, and 49.5 cm long), cast into 
trays, are positioned in the propellant section. The initial 
cross-sectional dimensions of the rectangular port (2.54 cm 
wide and 0.635 cm high), as shown in Fig. 3, are uniform 
along the entire motor. The ignition stimulus is provided 
by a fore-end, gaseous, pyrogen-type igniter, which discharges 
the ignition gases into the entrance section through a cen- 
trally located axial sonic nozzle. A rectangular replaceable 
nozzle, having a short convergent part and a well-defined 
throat, is attached to the motor aft-end. A detailed descrip- 
tion of the experimental facility is given in Chapter V. 

The first phase of the starting transient is the induc- 
tion interval (ignition lag) , which ends with the appearance 
of the first flame on the propellant surface. The delay 
associated with the processes .in the igniter depends upon the 
i articular ignition system and is therefore excluded from the 
general analysis. The beginning of the transient is deter- 
mined to be the onset of igniter discharge, defined as the 
start of pressure rise in the igniter chamber. A flow pattern 
with pressure, temperature, and velocity distributions devel- 
ops in the motor port at the igniter onset. Subsonic flow of 
the igniter gas with wall friction and heat transfer to the 
propellant and non-propellant port walls at low pressure (usu- 
ally between 1 and 3 atm for the low port-to-throat area 
ratios under investigation) is established in the chamber. 
Therefore, the motor nozzle may be either choked or unchoked. 
The heat flux to the propellant surface creates an unsteady 
temperature profile in the solid phase with constantly in- 
creasing surface temperature. 

The second phase of the transient is the flame spreading 
period which starts at first ignition and ends when the entire 
initially exposed propellant surface is ignited. The ignition 
criterion adopted in this study is that a point on the propel- 
lant surface ignites when it attains some critical ignition 
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(or autoignition) temperature, ' denoted here by , 

at which rapid runaway reactions leading to ignition^takft 
place. This ignition criterion is frequently used and leads 
tu results which are often indistinguishable from those ob- 
tained using other criteria. Many studies 2 - 4 5 >47 ,48 have 
shown that (under most conditions of motor ignition and for 
many propellants) Tp S/ ig to a good approximation is an in- 
variant quantity. It* should be noted, however, that for 
certain propellants it may be pressure-dependent. Thus the 
process of flame spreading along the propellant surface is 
described by the model of successive ignitions . 2 / 3 , 15 Once 
started, the flame spreading is accelerated by the increased 
heat flux, due to the rapidly increasing mass flow originat- 
ing from the already ignited propellant surface. The transi- 
tion from unchoked to choked flow in the port is considered 
in the analysis. Rapid chamber pressurization usually begins 
with the initiation of fast flame spreading. 

The third and last phase of the starting transient is 
the chamber filling period, which follows the completion of 
flame spreading. Severe erosive burning may take place dur- 
ing this period and thereafter, due to the high gas velocities 
in the chamber. Thus, the maximum press re, attained in the 
chamber, may be much higher than the desired equilibrium pres- 
sure. Significant pressure and temperature drops, as well as 
velocity increases, are established along the port. Follow- 
ing the pressure peak, a quasi-stea 'y situation develops and 
the head-end pressure decreases toward an equilibrium value, 
primarily due to diminishing erosive burning, as the port 
area increases. 

B. Basic Assumptions 

The following basic assumptions are used in the analyti- 
cal model: 

(1) All chemical reactions o<-cur on the propellant sur- 
face in a combustion zone, which is so thin that it 
can be considered as a plane. The combustion prod- 
ucts enter the main stream in the port with zero 
momentum. 

(2) The gases flowing in the port obey the perfect gas 
law. 

(3) Rate processes at the propellant surface are quasi- 
steady in the sense that their characteristic times 
are short, compared with the characteristic time of 
the pressure transient. 

(4) The propellant combustion products and the pyrogen 

igniter gas have the same values of c p , W and y. 
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(5) The flow in the chamber port is one-dimensional, 

with the fluid pressure, temperature, and velocity 
varying both with time and axial distance. 

Assumption (1) is based upon many theoretical and exper- 
imental studies on the combustion mechanisms of composite 
propellants. 41,95, 140 and is widely used by the investiga- 
tors . ^ / 95 , 141 A careful examination of this assumption is 
necessary, however, when this analysis is applied to double 
base and highly-metallized solid propellants. 

Assumption (2) is well acceptable for the usual range of 
high gas temperatures and moderate pressures in solid propel- 
lant rocket motors. 

Assumption (3) is justified for the following reasons. 
From Ref. 95 the characteristic time for surface reactions is 
T s ^ R u Tp S ap r /Ea r 2 . For typical values of composite propel- 
lant properties^ and burning rate varying from 0.3 to 0.6 
cm/sec at 20 atm, x s varies from 0.0005 to 0.002 sec, which 
is much lower than the smallest pressurization time observed 
in this study. 

Assumption (4) may not be needed, if the pyrogen igniter 
burns a solid propellant charge of the same composition as 
the propellant in the main motor. It has been proved by the 
solution of the analytical model that the effects of slightly 
different c p , W and y on the gas parameter distribution 
in the chamber are negligible during the induction interval, 
when they are most pronounced. The difference in temperature 
between the igniter gas and the propellant combustion gases 
is accounted for in the analysis, assuming instantaneous 
mixing with no chemical reactions. 

Assumption (5) is a good approximation to the real situ- 
ation in view of the highly turbulent flows encountered in 
HVT motors. Consideration of two-dimensional or three-dimen- 
sional phenomena, in addition to the variation with time, 
which exist to some extent in practical motor configurations, 
would make the solution of the problem very difficult to 
achieve with the existing mathematical tools. Willoughby^ 2 
has modified a computer program, recently published by Patan- 
kar and Spalding-*-^ f or a two-dimensional model of flow based 
on a turbulent compressible ooundary-layer theory, to calcu- 
late quasi-steady pressure drop along the solid propellant 
motor. He found a remarkable agreement between the pressure 
distribution calculated by Patankar's program and that calcu- 
lated by a simple one-dimensional model. 

Multi-dimensionality affects mainly particular processes 
such as the heat transfer to the surface, skin friction and 
erosive burning. Changes of flow properties across the bound- 
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ary layer are considered in the expressions for heat transfer 
and friction coefficient used in this analysis. The assump- 
tion of one-dimensional flow is consistent with the data re- 
duction methods used to obtain these correlations. 

C. Description of the Analytical Model 

A control volume element is established in the motor 
port, as shown in Fig. 4. The mathematical formulation of 
the approach to the above described problem consists of the 
following: 

(a) Mass, momentum, and energy conservation equations 
in unsteady, quasi-one-dimensional form for the gas 
phase; 

(b) Equation of state for the gas flowing in the motor; 

(c) Proper initial conditions at the start of the 
transient (onset of igniter flow) ; 

(d) Two boundary conditions at the fore-end of the pro- 
pellant section, obtained from a pair of ordinary 
differential equations, which describe the rate of 
change of pressure and temperature in the entrance 
section; 

(e) A third boundary condition, which describes the gas 
velocity at the entrance to the motor nozzle, for 
either choked, or unchoked flow; 

(f) Semi-empirical correlations for the convective heat- 
transfer and friction coefficients for the highly 
turbulent flow in the port; 

(g) Burning rate ’aw for the solid propellant, includ- 
ing the effects of initial temperature, pressure, 
and velocity (erosive burning) , and 

(h) A solid-phase heatup equation for determination of 
the propellant surface temperature during the in- 
duction interval, coupled to an ignition criterion 
for the solid propellant. 
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1. Conservation Equations for the Gas Phase 


The mass conservation equation for the control 
volume in Fig. 4, written for a unit volume, is 
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The momentum conservation equation, written for a unit 

control volume, is 
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The energy conservation equation, written in terms of 
the stored total (internal and kinetic) energy per unit mass, 
E, for a control volume of unit length, is 
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The unit-conversion factors g and J c are used in the 
conservation equations and throughout the analysis in order 
to carry out the entire study in a consistent, metric system 
of units, as denoted in the Nomenclature. 

The gravitational force is the only body force that ex- 
ists in the physical model analyzed in the study. However, 
the gravitational body force is negligibly small as compared 
to the surface forces for the physical situation considered, 
and therefore is not included in the mathematical model. Body 
forces such as centrifugal, axial acceleration, and electro- 
magnetic forces which are important in certain applications 
can be included in the model, when necessary. 

The above conservation equations, Eqs. (III-l) to (III- 
3), constitute a set of inhomogeneous, nonlinear partial 
differential equations. They are general in nature and can 
be applied to a broad range of practical rocket motors. 
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2. Simplification of the Conservation Equations 

For a practical solution of the conservation equa- 
tions a few simplifications are necessary. Each equation 
was nondimensionalized with respect to reference quantities 
specific for the general problem, and an order-of-magnitude 
analysis was performed to find the significance of each term. 
This procedure is described in Appendix A. As a result of 
the order-of-magnitude analysis the effects of the following 
terms were neglected: 

(1) Rate of mass accumulation in the free volume, 
created by the propellant surface regression [as 
compared to the rate of mass addition due to burn- 
ing, last term in Eq. (III-l)]. For the extreme 
operational conditions of p = 100 atm, W = 40 
g/g-mole and T = 1500°K, there is p = 0.02 Pp r 
(if pp r = 1.6 g/cm^) . For the usual case of 

p = 40 atm, W * 22 g/g-mole and T = 2500°K, 

p = 0.0026 pp r ; 

(2) Forces between gas molecules, due to the viscous 
stress in the axial direction, created by the axial 
velocity gradient [2nd term on the right-hand side 
Of Eq. ( I X I — 2 ) ] ; 

(3) Axial heat conduction between gas molecules [1st 
term on the right-hand side of Eq. (Ill— 3) } ; and 

(4) Rate of work done by the internal viscous forces 
and viscous dissipation [3rd term on the right-hand 
side of Eq. (III-3)]. 

The order-of-magnitude analysis shows that the friction 
force term in Eq. (III-2) cannot be neglected in laboratory- 
size HVT motors, even though it is small. Indeed, the para- 
metric study conducted to investigate the effects of differ- 
ent parameters has confirmed this conclusion, as described 
in Chapter VI. Another reason for retaining the friction- 
force term is the tendency to have an analysis, which is able 
to check the effect of pre-firing propellant-surface roughen- 
ing. This roughening has been explored as a means of short- 
ening ignition delay times. However, it does not lead to 
reproducable results and is not widely used. 

Many practical rocket motors are designed with initial 
port area which varies with distance from the fore-end. It 
should be pointed out that the analysis considers this ini- 
tially prescribed function Ap f ^(x), as well as the varia- 
tion of port area after ignition Ap(x,t), caused by the 
uneven axial distribution of burning rate along the motor 
port. It should be also mentioned that the analytical model 
corresponds closely to motors with pyrogen igniters, but can 
be readily modified for application to motors using pyrotechnic 
igniters . 
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After deleting the above-mentioned terms of secondary 
importance, the conservation equations are rewritten in the 
following. 

The mass conservation equation is 
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Following the substitution of the mass conservation 
equation, Eq. (III-l) , into Eq. (III-2) , and the use of a 
conventional definition for the friction coefficient (f = 
2gx w /pu 2 ) , the momentum conservation equation becomes 
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After the substitution of the mass and momentum conser- 
vation equations [Eqs. (III-l) and (III-2)] into Eq. (III-3) 
(performed before the nondimensionalization and order-of- 
magnitude analysis) and some algebraic manipulations, the 
energy conservation equation, in terms of the static gas 
temperature, becomes 
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For simplicity, the specific heat at constant pressure, 
Cp, the mean molecular weight, W, and the ratio of speci- 
fic heats, y, are assumed to be constant throughout the 
entire starting transient. Also, species diffusion effects 
due to concentration gradients caused by mixing of the pro- 
pellant combustion gases and the igniter gas, and by axial 
variation of the flowing gas parameters, are neglected. The 
pressure-dependence of the propellant adiabatic flame temper- 
ature is neglected too. 

3. The Governing Equations and Their Initial Conditions 

After substitution of the perfect gas law, p = pRT, 
into Eq. (IJI-5) , the momentum equation reduces to the first 
governing equation: 
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Two other governing equations are obtained by substitu- 
tion of the perfect gas law into the mass and energy conser- 
vation equations [Eqs. (III-4) and (III-6) ] , combining the 
two, and some algebraic manipulations: 
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The three governing equations [Eqs. (III-7) , (III-8) , 
and (III-9) ] are first order, inhomogeneous, coupled, and 
nonlinear partial differential equations. The independent 
variables are time t, and distance along the port x, and 
the unknowns are the gas velocity u, temperature T, and 
pressure p. This set of governing equations contains all 
essential gas-dynamic features of the physical model, as 
applied to propellants with no :>r low aluminum content. How- 
ever, rocket motors that use highly aluminized propellants 
that do not produce large AI2O3/AI agglomerates (i.e., the 
aluminum burns very close to propellant surface) are also 
closely approximated by this set of equations. When the 
aluminum particles are sufficiently small, the two phase flow 
effects are not important and effective values of y, cp, and 
W enable the aluminum/combustion gas mixture to be consid- 
ered as if it were gas only. 


The three necessary initial conditions, which describe 
a state of uniform distribution in the propellant and aft 
sections of the motor, are: 


u (0,x) = 


T ( 0 ,x) 


T i 
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where uj., T^, and p^ are the initial gas velocity, tem- 
perature, and pressure, respectively, in the motor. To facil- 
itate the start of the numerical solution, it is assumed that 
Pi is slightly higher than the ambient pressure (by 0.0001 
atm) , with resulting small velocity ui in the port. 

Examination of the radiative heat transfer from the 
cloud of igniter gases or propellant combustion gases in the 
chamber to the port walls-^4 ,145 shows that for the case and 
conditions considered in this study with unaluminized propel- 
lant, heat transfer by radiation is less than 5% of the con- 
vective heat transfer. This is in agreement with the measure- 
ments and findings of many investigators. 7 ' 56 /1Z0 

The following expression for the rate of heat transfer 
to the propellant and nonpropellant walls, q^, pertinent 
to the physical and experimental configuration (see Fig. 3) 
is used in the analysis: 
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where h C p, and h pw are the local convective heat -transfer 
coefficients over the propellant and nonpropellant port walls, 
respectively. 

In particular cases such as the presence of solid parti- 
cles in the igniter gases and the use of metallized propel- 
lants, the radiative heat transfer may become substantial 
and should be considered. 4 


The Physical Boundary Conditions 


The physical boundary conditions are specified in 
accordance with the flow conditions at both ends of the motor. 


A comprehensive study of the structure (boundary and 
internal shock-pattern) of axisymmetric free jets exhausting 
into still air and into supersonic streams was conducted by 
Love and co-workers. 1 ^ A research group in Rocketdyne in- 
vestigated experimentally the gas-dynamics and heat-transfer 
characteristics of simulated head-end igniters in configura- 
tions generally identical with that in this study. 50# 52, 53 
Figure 5, taken from Ref. 53, shows the jet structure of a 
head-end sonic igniter, exhausting into a duct with a nozzle 
at its other end, for two possible cases: (1) expansion into 

a duct with relatively small diameter (see Fig. 5a) , and (2) 
expansion into a duct with relatively large diameter (see Fig. 
5b) . The sizes of igniter nozzle and motor port used in this 
study imply the existence of the case, described in Fig. 5a, 
in vertical direction, whereas the case, described in Fig. 5b 
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is valid in the horizontal direction (see Fig. 3) . Therefore, 
the three-dimensional flow pattern in the entrance section 
of the experimental motor during the igniter operation is 
quite complex. 

In the analytical model of this study it is assumed that 
fast stirring takes place in the entrance section, and there- 
fore the gas properties are unifoim therein. The mass and 
energy conservation equations for a control volume, taken to 
be the entire volume of the entrance section, provide two 
boundary conditions at the fore-end of the propellant section 
as described in the following paragraphs. 

The mass conservation equation for the gas in the entrance 
section is: 


V 


dp es (t) 


es 


dt 


* V * 1 ' A p,« P es< t ’V U 'Bs (t 'V (III - 12 > 


where mig(t) is the igniter mass flow rate into the entrance 
section. y The second term on the right-hand side represents 
the mass flow rate out of the section. 


The energy conservation equation for the entrance sec- 
tion is 
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Order-of-magnitude evaluation, confirmed by the results 
of the numerical solution, shows that the rate of change of 
the kinetic energy in the entrance section is negligible in 
comparison with the rate of change of the thermal energy. 
Hence, the second term in the parenthesis on the left-hand 
side of Eq. (III-13) can be ignored in this analysis. 


After substitution of the equation of state into Eq. 
(III-13) , and rearrangement, the following equation is obtained 
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Substituting the equation of state and Eq. (III-14) into 
Eq. (III-12) one gets the following equation after some alge- 
braic manipulations: 
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The igniter mass flow rate, mig, is a prescribed func- 
tion of time, determined from the experimental study, as 
described in Chapter V. High heat-transfer and shock-pattern 
losses, taking place in the entrance section, are found^'53 
to lower the effective mean temperature of the igniter gas 
there, T ig/ which is assumed to be constant during the 
transient. This point is discussed in Section D of this 
Chapter . 


Equations (III-14) and (III-15) form a pair of coupled 
ordinary differential equations, which describe the rate of 
change of the gas pressure and temperature in the entrance 
section. Their solution provides two boundary conditions for 
the governing equations, Eqs. (III-7) to (III-9) . The two 
initial conditions for this pair of equations are 

p es <0) " p es,i and T es (0) * T e s,i UII-U) 


The physical boundary condition at the aft end of the 
motor chamber is derived assuming quasi-steady isentropic 
flow between the entrance to the motor nozzle, which coin- 
cides with the chamber end, and the nozzle throat (see Fig. 
2). This assumption is justified by the very small length 
of the nozzle entrance. Using mass and energy conservation 
equations (as well as isentropic relations) between the two 
sections, the following equation is derived: 
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When the motor nozzie (with no divergent section) is 
unchoked (which is the usual case during the induction inter- 
val for very low port-to-throat area ratio and low to moder- 
ate igniter mass flow rates) , the pressure at the nozzle 
throat, p t , equals the ambient pressure, p a , and Eq. (III- 
17) becomes 
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It should be noted that Eq. (III-18) is given incorrect- 
ly by Shapiro on p. 966 of Ref. 148. This is probably a 
typographic error, since the plots accompanying the formula 
(p. 965 of Ref. 148) correspond to the correct equation, 

Eq. (III-18) . 

For the case of a choked motor nozzle, Eq. (III-17) be- 
comes the following implicit relation: 
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The analysis considers both the unchoked and choked 
flow cases. The computing scheme in the numerical solution 
checks constantly the ratio of ambient pressure to stagnation 
gas pressure at the nozzle entrance for the choking condition 
and either Eq. (III-18) or Eq. (III-19) is used, accordingly. 
Usually the nozzle is choked at the early stages of the flame 
spreading phase. For instance, for m^g = 15.0 g/sec at 
A p /A t 11 1.2 the motor nozzle became choked when 10% of the 
exposed propellant surface had ignited. 

It is noteworthy to mention, that if one considers a 
case in which a diaphragm is placed at the nozzle entrance, 
or a plug is inserted in the nozzle throat, the aft-end 
boundary condition for the period of time preceding the dia- 
phragm burst, or the plug blowout, is simply 


u E (t) « 0 


(III-20) 
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It should be also noted that a careful consideration of 
the physical situation after the motor nozzle throat is need- 
ed. Existence of a diverged part and an exhaust system 
should be taken into accour. . when establishing the aft-end 
boundary condition. 

In Appendix B the governing equations are presented in 
a form which shows the effect of mass addition, unsteadiness, 
heat transfer, etc. on the axial variation of pressure, tem- 
perature, and velocity. That presentation of the governing 
equations clarifies the contributions of this study as com- 
pared with the quasi-steady one-dimensional analyses. 

D. Correlations for uhe Convective Heat-Transfer Coefficient 


The flow in the port of HVT motors is highly turbulent 
(typical Reynolds numbers are given in Table 4) . Dittus- 
Boelter's correlation for the convective heat-transfer co- 
efficient for a fully-developed turbulent flow in a smooth 
pipe is widely accepted for use in so lid -propellant rocket 
motors. This correlation has the following form: 


Nu d = 0.023 Re^’Sr 0,4 


(III-21) 


where subscript d denotes quantities based on duct diameter, 

For turbulent flow in ducts with non-circular cross- 
section (and for Pr > 0.6 and Re^ > 7,000) the diameter 
successfully used in Eq. (III-21), as in expressions for the 
friction coefficient, is the equivalent hydraulic diameter 
dh*^ 2 3 4 ^'^ 0 The applicability of the concept of the hydraul- 
ic diameter for turbulent flow has been verified experimen- 
tally up to a Mach number of 1.0.150,151 

The flow in a solid-propellant rocket motor with an in- 
ternal-burning grain differs from a fully-developed turbulent 
flow in a smooth duct. Therefore, if a proper application of 
Dittus-Boelter ' s correlation is to be made to a practical 
rocket motor, corrections for the following effects and phen- 
omena must be considered: 

(1) Entrance effects at the fore-part of the motor, 
resulting from the distance necessary for the de- 
velopment of the hydrodynamic and thermal boundary 
layers ; 

(2) Sharp variation of the fluid properties across the 
boundary layer; 

(3) Highly transient I low conditions in the port; 

(4) Propellant surface roughness which increases the 
heat transfer to the surface; 











(5) Decomposition-gasification of the unburned propel- 
lant surface, which attenuates the heat transfer 
process ; 

(6) The general configuration of a solid rocket motor 
with head-end ignition, in which a sonic jet is 
discharged from an orifice of a relatively small 
diameter into a duct that is constricted by an aft- 
end nozzle [in contrast to the configuration for 
which Eq. (III-21) applies directly] ; 

(7) The configuration and stagnation properties of 
the igniter jet, or jets, and the motor pressure 
distribution. 

(I?) Impingement of incandescent particles, which may 
exist in the flowing gas, on the unignited propel- 
lant surface; 

(9) Changes in flow conditions downstream of the igni- 
tion front during the flame spreading, and 

(10) Heat-releasing or heat-absorbing processes on the 
unignited solid-propellant surface, such as gas- 
phase condensation, and exothermic and/or endotherm- 
ic reactions, which combine with the heat-transfer 
from the hot gas flow. 

Consideration of all these effects is a very complicated 
task, especially because some of them are unclear quantita- 
tively. Therefore, an experimental determination of the 
heat-transfer correlation is necessary, even though it is 
difficult to simulate all practical conditions for this pur- 
pose. 


A brief discussion of the above-mentioned phenomena af- 
fecting the heat-transfer coefficient is pertinent here, be- 
cause of its importance in calculating the starting transient 
of HVT motors. 

The entrance effects have been observed by many investi- 
gators . 15 ,47 , 52 , 53 , 152-154 The developing flow at the en- 
trance part of a duct, or a motor, has been correlated by 
several methods: (a) assuming laminar or turbulent flow, 

according to the flow conditions, over a flat plate^ >56 , 61 ,77 
(thus the coefficient in Eq. (III-21) is increased and the 
effect of length appears through basing Reynolds and Nusselt 
numbers on distance) ; (b) multiplying the Dittus-Boelter 

expression by a power function of the length-to-diameter 
ratio, such as ai(x/d) s , where ai takes values between 1 
and 2, and the values of s range between -0.05 and -0.2; 
155-157 anc j ( c ) multiplying Eq. (III-21) by an exponential 
function of the length-to-diameter ratio, usually of the 
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form & 2 [1 + a 3 exp (- a 4 X r /d) ] , ' where a^, a 3 , and 

ate constants. All three methods result in similar correc- 
tions, with pronounced entrance effects up to 6 to 20 duct 
diameters downstream. Maximum heat transfer for pyrogen 
igniters fired into instrumented ducts was found to occur in 
the igniter jet impingement zone and is up to 6 times that 
calculated by Eq. (III-21) . 

A consideration of the variation of fluid properties 
between the bulk of flow and the propellant surface is 
necessary in rocket motors because of the large temperature 
difference and considerable cross-sectional variation of 
transport properties. Several empirical methods of correc- 
tions to Eq. (III-21) have been proposed for this purpose: 

(a) multiplying Dittus-Boelter ' s correlation by a power 
function of the gas-to-surface temperature ratio of the type 
a 5 (T, ^g) 3 , or of the ratio of gas viscosity at the bulk turn 
perature to gas viscosity at the surface temperature; 2 
155,159,1 (b) evaluating the physical fluid properties in 

Eq. (III-21) at an average film temperature T a f « (T + T s )/2 
and using the bulk velocity, as correlated successfully by 
Humble, et al 154 for subsonic flow in long tubes and by 
Bartz^ 61 for highly turbulent flow in convergent-diverge-.c 
nozzles, and as recommended by McAdams. 1 ^ 

The effect of flow unsteadiness on 'he convective heat 
transfer is. in general, not well known. ' 63 As described 
by DeSoto and Friedman, »t onset of flow over a body, 
the boundary layer is not distinguishable and the entire 
flow pattern appears inviscid. As time progresses a boundary 
layer materializes. Thus the convective heat-transfer co- 
efficients should be higher during the transient period than 
under steady-state conditions. Carlson and Seader^®'^ 

(using hot nitrogen at pressures from 570 to 620 psig and 
temperatures from 550 to 950°F) found a strong time dependence 
of the heat-transfer coefficient h c . In some tests, h c 
decreased by a factor of two at the end of 0.8 second at the 
location of the beat transfer maxima. In a later report, 
Wrubel and Carlson 53 attributed the transient effects observed 
by Carlson and Seader to undetected conduction losses and they 
considered time dependent effects to be negligible. Experi- 
ments 53 with high driving temperature potential (by firing 
pyrogen igniters into an instrumented duct) did not show 
transient heat-transfer behavior, except in areas with appre- 
ciable axial conduction losses. Knuth^ 62 showed from a 
theoretical analysis that for most applications forced-con- 
vection heat transfer rates to surfaces with time-depe dent 
temperatures can be calculated to a good approximation using 
equations developed for the steady-state case. A comprehen- 
sive experimental worlj conducted at the University of Utah 
and reported in 1965 5- ' 64 has confirmed the conclusions of 
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Knuth. In a later study, Baer and Ryan observed unsteadi- 
ness effects (see Chapter II, Section A). Sparrow and 
Siegel^” treated analytically the transient nature of the 
convective heat-transfer coefficient. 

The effect of piopellant surface roughness in increas- 
ing the convective heat transfer and thus decreasing ignition 
delay times and increasing flame-spreading rates has been 
widely observed. Unfortunately, well established quantita- 
tive data is missing. The increasing amount of experimental 
data indicates that the heat-transfer coefficient increases 
in a much slower rate with increasing surface roughness, than 
the friction coefficient. ^ >159, 164-166 Recent studies!67-170 
have developed semiempirical correlations for heat transfer 
to rough surfaces in subsonic and supersonic turbulent flow 
and have concluded that the Reynolds analogy is invalid for 
rouqh surfaces. This invalidity increases with Reynolds num- 
ber .1^6 , 169 

The attenuation of convective heat-transfer coefficient 
and skin friction for turbulent flow over a surface with 
blowing has been the subject of many theoretical and experi- 
mental investigations.^-!' >118, 171-179 However, until the 
physical law of the decomposition-gasification process taking 
place on the propellant surface is well established, appli- 
cations to solid rocket motors can be done only in a specula- 
tory manner. 

The uniqueness of a soiid-propellant rocket motor as 
compared to a pipe, as far as the convective heat transfer 
is concerned, was shown clearly by the studies of Carlson 
and Seader-0 - -2 and Wrubel and Carlson, 53 with the addition- 
al interpretation of Cohen, et al.1^8 These studies have 
revealed the effects of igniter-nozzle configuration, motor- 
port to igniter-nozzle diameter ratio, and motor port-to- 
throat area ratio on the heat transfer distribution for 
relatively short motors. The importance of the igniter jet 
structure has also been shown. The secondary effects caused 
by the port-to-throat area ratio are primarily related to 
igniter jet expansion characteristics. It appears t.iat these 
parameters have to be considered in HVT motors, since the 
above-mentioned studies were conducted with low port-to-throat 
area ratios (1.1, 2.0, and 3.0). This conclusion has been 
supported by the results of this investigation. 

5 

Most and Summerfieid studied the effect of molten alum- 
inum particles liberated from the ignited portion of the grain 
during the flame spreading on the ignition of the remainder 
of the surface. No significant effect was observed. The 
parti ^les that hit the unignited surface neither attached 
themselves to the surface, nor caused any spread of ignition 
as they traveled downstream. It was concluded that the quan- 
tity of hot aluminum particles in the flow was not sufficient 
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to augment significantly the flame spreading. However, radi- 
ation from densed clouds of incandescent particles should be 
considered. 

4 5 

The experimental study of Jensen and Cose did not find 
any evidence that postulated vortex flow ahead of the advanc- 
ing flame affects significantly the heat transfer downstream 
of the flame front. Mostly also corcluded from his heat-flux 
measurements that the convective heat transfer during flame 
spreading is only very weakly dependent on the relative loca- 
tion of the flame and he related it to the original leading 
edge. However, when high burning rate propellants are used, 
with resulting large blowing parameter, this point should be 
considered carefully. 

Propellant surface reactions have been considered in 
some heat-transfer studies* 4 ' 45 whereas very little has been 
published on gaseous condensation phenomena. 2 

After reviewing the large amount of literature on the 
convective heat transfer in real and simulated solid-propel- 
lant rocket motor, and in view of the highly turbulent flow 
in HVT motors in general, a correlation for the heat-transfer 
coefficient was deduced in this study. This deduction was 
mac. !, keeping in mind the features of the experimental motor, 
for which the analytical model was first solved. 

The expression for the local heat-transfer coefficient 
is based on the Dittus-Boelter correlation, Eq. (III-21) , 
the dimensionless parameters, in which, are calculated using 
the local flow properties. This coupling to the chamber gas 
dynamics is necessary due to the large velocity, temperature, 
and pressure variations along the port, which are calculated 
by solving the governing equations, Eqs. (III-7 to III-9) 

(see Chapter VI) . Corrections to Eq. (III-21) are introduced 
for the entrance effects, cross-sectional variation of gas 
temperature, and port-to-throat area ratio. 

The corrections for the effects of length-to-diameter 
ratio and the variation of gas temperature from the core of 
flow to the surface are based on the studies of Humble, et 
al,* 54 cholette,* 55 and Bartz.*°l The correction for port- 
to-throat area ratio is base' on the interpretation of Cohejj, 
et al* 5b to the studies of Carlson, Seader, and Wrubel. 5 0-5- 
The latter contain the best available experimental heat- 
transfer data concerning HVT motors with head-end pyrogen 
igniters . 
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Introducing the above-mentioned corrections 
the following correlation is obtained: 
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where : ag is a constant; x^- is the distance from the 

point of reattachment of the igniter jet, and the physical 
properties of the gas (except the velocity) in the Nusselt, 
Reynolds and Prandtl numbers are evaluated at the average 
film temperature T a f. 

From Eq. (III-22) an expression for the local heat- 
transfer coefficient developed for this study is obtained in 
the following form: 


h c = 
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(x r d h r 0 ' 1 ( V A t )0 


(III-23) 


The temperature dependence of the viscosity is deter- 
mined from the recent calculations of Svehla-^O for the vis- 
cosity of air at high temperatures. It was shown by Bartzl^l 
that Svehla's values are very close to those in the previous- 
ly published NBS data. 1^2 

Over the range from 1000 to 3500°K the calculated high- 
temperature viscosity values for airl^O were fit in this 
study by the formula 

y = 0.87 x 10' 6 W 0,5 t(°K) 0,65 poise (III-24) 


Prandtl number is assumed to be constant and calculated 

from Svehla's equation^O 


Pr 


Y _ 
1.77y-0. *5 


( III-25 ) 


which gives approximately the same value of 0.71 for y = 1.22 
to 1.26 (usual propellant combustion gases). 

Expressing the gas density p a f by the local pressure 
and average film temperature 


P af = p/RT af (III-26) 

and substituting Eqs. (III-24) and (III-26) into Eq. (III-23) 
the following expression is obtained: 


. 0 • 6„0 • 1 , / n t 0 . 8 m —0 .67 . —0 .1,, /, .0.4 

(III-27) 
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where a^ = (0.87 x 10 ^)®*^ag. 
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As mentioned already, the best available pertinent ex- 
perimental data on the convective heat transfer are those 
obtained by Carlson, Seader and Wrubel 5 ^“53 at Rocketdvne. 
Examining carefully these data, Cohen, Derr and Price^*® 
proposed tne following correlations for the convective heat 
transfer from an axial igniter with a sonic nozzle: 

At the region of igniter- jet impingement 
(maximum heat transfer) : 

Nu d = 0.0615 Re d 0 ' 8 (A /A fc ) 0,4 (III-28) 

Upstream of the impingement point: 

Nu, = 0.7 Re, 8 ' 8 (A /A, ) 8 ' 4 (0 . 022 + 0.066 x/x. ) 
a a v p t lm 

(III-29) 

Downstream of the impingement point : 

Nu, = 0.7 Re, °‘ 8 (A /A. ) 8 ' 4 [0 . 031 + 0.057 exp(-0.5 x /d) ] 
d d p t r 

(III-30) 


Unfortunately, the data in Refs. 50 through 53 extends 
only to x/d =10 approximately. Eq. (III-30) gives almost a 
constant value for Nu<j after that point. The experimental 
motor in this study extends up to x/d^ = 62 with the propel- 
lant leading edge being at Xp/d^ = 6.06 (or Xp/dj^ = 5.06 
on the average) . The constant ag in Eq. (III-23) was esti- 
mated so that over the propellant lead! lg edge and up to 
x/dh = 10 the heat-transfer coefficient calculated by Eq. 
(III-27) will be equal to that found in Refs. 50 through 53 
and calculated by Eq. (III-30) . A value of ag = 2.56 x 10“2 
was determined. Therefore, the expression for the lor'l heat 
transfer coefficient [Eq. (III-23)] developed in thi study 
becomes 


h c = 1.56 x 10“ 3 c p Pr" 8 * 6 W 8 ' 1 (pu/R) °‘ 8 T af ° ' 67 ( x r d h ) “° * 1 < A p /A t ) 8 * 4 

(III-31) 


The position of igniter jet impingement, x^ m , is found 
by the jet-boundary calculations presented by Love, et al. 47 
For most of the experiments, the position was in the entrance 
section approximately at one port diameter downstream of the 


igniter nozzle exit, i.e., 


x im & 


d h . 


Expression (III-31) was used throughout this study. It 
is based on experimental studies of others and it is recommend- 
ed that experimental investigation will be conducted in the 
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test motor for verification (see Chapter VII) . Good agree- 
ment between measured and calculated pressure-time-distance 
curves was obtained using Eq. (III-31) . Since the computer 
program is flexible, there is no problem in using other 
expression, proved to be more precise, in the numerical solu- 
tion of the analytical model. 

Before local ignition, both the local heat-transfer co- 
efficients to the propellant and nonpropellant walls of the 
port, h cp and h cw respectively, are calculated by expres- 
sion (III-31) , i.e., h C p = h cw ~ h c . After local ignition 
h C p = 0 and h cw = h c • 

Equations (III-28 to III-30) are used to calculate the 
effective mean temperature of the igniter gas, Tj_. Calcu- 
lation of the stagnation temperature in the igniter combus- 
tion chamber using pressure measurements during steady state 
igniter operation shows 5% losses due to incomplete combus- 
tion and heat transfer, as compared to the adiabatic flame 
temperature. The high heat transfer in the entrance section 
causes an additional drop of between 25 and 35%, according 
tc the port-to-throat area ratio. It should be mentioned 
that Most, 15 using similar ignition system with convergent- 
divergent igniter nozzle and an entrance section a little 
shorter than the one in this study, measured igniter gas 
temperatures over the propellant leading edge, which are 
about one half of the igniter adiabatic flame temperature. 

E. Correlation for the Friction Coefficient 


The correlation for the friction coefficient, f, is 
deduced from Colebrook's expression for turbulent flow in 
pipes with roughness , 150 , 183 ,184 v/hich is the basis of 
Moody diagram. Colebrook's formula is 


1 _ 

/? 


- 0.8686 


In 


e c /d, 
s' h 

3.7 


2.51 

Re,/? 

d 


(III-32) 


where e s/^h I s the relative equivalent sand roughness. 

Observed effects of distance from the entrance^^ »185 
and cross-sectional variation of fluid properties are taken 
into account. Humble, et al^54 correlated successfully 
friction-coefficient measurements in smooth tubes with 
Prandtl ' s universal law of friction for smooth pipes evalu- 
ating density and viscosity at the average film temperature. 
They did not correct for observed entrance effects. Shapiro 
and Smith 1 ° 5 observed entrance effects up to 20 to 30 tube 
diameters from the tube inlet. Therefore, Reynolds number 
in Eq. (III-32) is calculated with the local flow velocity, 
whereas density and viscosity are evaluated at the local 
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\ 


i ' 

i 


0 

0 



pressure and average film temperature. Entrance effect is [ 

considered by multiplying the friction coefficient by a power [ 

function of the distance-to-diameter ratio as for the heat- i 

transfer coefficient, a^ (x r /d^) ^ • 1 , and evaluating the coef- f 

ficient so that this factor will become unity at twenty duct [ 

diameters. 1 

After introducing the aforementioned corrections into 
Eq. (III-32) , and some algebraic manipulations, the implicit 
correlation for the friction coefficient developed in this 
study is: 


f 



/ d h 

77 


0.449(d h /x r ) 0,1 

1.27 RW°' 5 T af 1,65 (d h /x r ) 

— ^ p ud ^ f A . 5 


0.05 


(III-33) 


Friction is considered over the entire surface during 
the induction interval and only downstream of the ignition 
front after first ignition. Zero wall friction is assumed 
at the burning surface due to the large friction attentuation 
caused by the blowing. Since the inert perimeter is small 
as compared to the burning perimeter, it is assumed in the 
model for simplicity that after ignition at a position the 
local friction coefficient is zero, i.e., the friction over 
the entire perimeter is neglected. 

F. The Burning Rate Law 

The solid propellant used in this study is a composite 
PBAA-AP propellant, composed of 20% PBAA-EPON 828 and 80% 
bimodal ammonium perchlorate. The properties of this pro- 
pellant are listed in Table 1. The non-erosive burning rate 
law, of the Saint-Robert 1 s type, as determined by strand- 
burner measurements is 

r Q (cm/sec) = 0.927 x 10 _2 p (gf/cm 2 ) 0,4 (III-34) 

Figure 6 shows this burning rate law in the onventional 
log tq vs log p form. Pressure is given in gf, cm 2 for 
consistency of units in the entire analysis. It is recalled 
that 


10^gf/cm 2 = 0.968 atm * 0.981 bars = 14.22 psi 

Large enhancement of the burning rate due to erosive 
burning, encountered in HVT motors, is taken into account by 
adopting the widely-used and accepted Lenoir-Robillard 1 s semi- 
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empirical burning rate law 112 

r = ap n + kh exp(-8rp /up) (III-35) 

where the non-erosive component is given by Eq. (III-34) , the 
local zero-blowing heat-transfer coefficient, h c , is calcu- 
lated according to Eq. (III-31) , and k and 8 are constants. 

This is actually a modified form of Lenoir-Robillard * s 
law, since the expression originally proposed for h c is 
based on Chilton-Colburn ' s correlation for turbulent flow over 
a flat plate, whereas the expression used here [Eq. (III-31) ] 
is a corrected Dittus-Boelter * s correlation. The local burn- 
ing rate as given by Eq. (III-35) is thus strongly coupled to 
the chamber gas dynamics and the propellant surface temperature. 

The erosive-burning constant, k, is estimated from 
semi-theoretical considerations as proposed originally by 
Lenoir and Robillard.l 1 ^ The erosive-burning exponent, 8, 
is evaluated from burned distance versus axial position 
measurements from propellant recovered from water-quench ex- 
periments conducted in this study. It was also adjusted to 
fit the experimental pressure-time curves for all values of 
Ap/At tested. The values of k and 8 used in this inves- 
tigation are listed in Table 1. 

The high pressurization rates encountered in HVT motors 
suggest the consideration of dynamic burning effects. These 
effects may be of two kinds: (1) due to burning of a pre- 

heated propellant layer on ignition, a phenomenon frequently 
oversized by undersized igniters, and (2) due to burning of 
a propellant layer, the temperature profile in which is not 
in phase with the rapidly varying pressure in the gas phase. 
Dynamic burning is not considered in this study. Its inclu- 
sion into the current model would further improve the analy- 
sis and the agreement with the experimental results. Dynamic 
burning is further discussed in the outline of suggestions 
for future work (see Chapter VII) and in Appendix C. 

G. Determination of the Propellant Surface Temperature 

To determine the amount of heat transferred to the solid 
propellant in the unignited region and the flame front loca- 
tion according to the ignition criterion, the propellant sur- 
face temperature must be evaluated as a function of time and 
distance along the port. The one-dimensional heat equation 
for the solid phase at a fixed axial location is 

3T (t,y) 3 2 T (t,y) 

TT = OL 

9t P r 3y 2 


(III-36) 


where y is a coordinate normal to the surface with its 
positive direction into the solid (see Fig. 7) . The propel- 
lant is assumed to be inert until the critical ignition tem- 
perature is reached. 


The initial condition for Eq. (III-36) is 

T (0,y) = T . 
pr v pi 


(III-37) 


and 


The boundary conditions for the heat equation are: 

Tpr (t,“) = T pi (III-38) 


9T h (t) 

-JF<t,0> = - — y [T ( t) - T It)] (III-39) 

ay a ps 


The second boundary condition specifies that at the pro- 
pellant surface the rate of heat transfer from the gas phase 
to the solid is equal to the rate of heat conduction into the 
solid, with the variables h c (t) and T(t) depending on the 
local flow parameters. It follows that Eq. (III-39) is cou- 
pled to the gas phase conditions. This implies that the heat 
equation must be solved simultaneously with the governing 
equations for the gas phase to yield the time -dependent tem- 
perature profile within the solid at each given axial posi- 
tion in the motor. However, the main interest in this analy- 
sis is not a detailed knowledge of the temperature profile 
inside the propellant slab, but rather the propellant surface 
temperature, which determines the rate of heat transfer to 
the unburned propellant and the flame spreading rate. 

A conventional way to treat the parabolic partial differ- 
ential equation (III-36) is to reduce it to an ordinary dif- 
ferential equation by Laplace transformation and solve for 
the surface temperature. 186 However, since no explicit ana- 
lytical expressions for h c (t) and T(t) are available, 
the analytical solution can only be stated in the integral 
form. Solving for the propellant surface temperature along 
the port at each time step requires a lot of computations 
and computer storage spaces. 

An approximate solution of the heat equation, which is 
quite accurate and saves much computation time may be obtain- 
ed by using the so-called integral method. 18 '' 188 Applying 
this method, a quantity 6 (t) (called the penetration dis- 
tance) is defined, so that for y 1 6 the solid is at equil- 
ibrium temperature and there is no heat transfer (see Fig. 7) . 
It is assumed that a third-degree polynomial in y represents 
the transient temperature profile in the solid: 
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T pr = k Q (t) + k x (t)y + k 2 (t)y 2 + k 3 (t)y 3 (HI-40) 

where the coefficients kQ, k^, k 2 and k 3 may be time- 
dependent. To determine them, four boundary conditions are 
necessary. These are: 

(1) T pr (t ' 6) = T pi (III-41) 


( 2 ) 


3T 


-^■(t,6) = 0 (no heat transfer at y = p) 


9y 


(III-42) 


( 3 ) 


2 

3 T 


(t,6) = 0 (the smoothing condition) 

3y^ (III-43) 


(4) 


9T pr h c 

= - jS_!T<t> - T ps (t)] 


(III-39) 


Solving for the coefficients krj to k 3 and substitu- 
ting into Eq. (III-40), the following expression for the 
temperature profile is obtained: 

3 


T = T . + 
pr pi 


h c (T - T ps )(6 " y) 

£ ~r“ 


3A 

pr 


(III-44) 


= T . + (T -T . ) 
pi ps pi 


<$-y 


!ni-',5) 


Figure 7 shows the temperature profile in the solid 
expressed by Eq. (III-45) . 

The surface temperature is obtained by substituting 
y = 0 into Eq. (III-44) 




v T -y 6 


3A 


pr 


(III-46) 


Following the integral method, a heat balance integral 
is written for the penetration distance, 6(t): 
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0 ± j 

l pr p*- = at < T pr- T pi )dy 

0 0 


(III-47) 


Equation (III-47) is analogous to the momentum integral 
in boundary layer theory. The left-hand side becomes, in 
view of conditions (III-39) and (III-42) , 


- v 


a h 


^,T-T ps , 


(111-48) 


The right-hand side of Eq. (III-47) becomes, using 
Eqs. (III-44) and (III-46) , 


ar (V T pi ,dy = 


3, a VV 

T*pr 3£ K^T-T ps ) 


{III-49) 


The rates of change of the heat-transfer coefficient, 
h c , and gas temperature, T, with respect to time, are 
neglected as compared with the rate of change of T ps , which 
is completely justified for the entire induction interval, 
and the major part of the flame spreading period. Performing 
the differentiation in Eq. (III-49) and equating the result 
to Eq. (III-48) the following ordinary differential equation 
for the propellant surface temperature is obtained: 


dT ps . 4 yh c 2 (T-Tp.) 3 

dt 3 V 2<T ps- T pi )(2T - T ps-V 


(111-50) 


To avoid the singularity at t = 0, the initial condi- 
tion for Eq. (III-50) is taken as 


T ps (0 > " T pi + c 


(III-51) 


where e has a small value, e.g., 0.1°K. 
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18 8 

It was shown by Goodman, that for constant external 
heat flux the surface temperature increase of a semi-infinite 
slab calculated by the integral method using a cubic tempera- 
ture profile is only 2% higher than that calculated from the 
exact solution of Eq. (III-36) . As a check of Eq. (III-50) , 
the heat transfer coefficient (h c ) was varied linearly with 
time in a manner that approximates the actual heating rates 
in the physical model. Solution of Eq. (III-50) by a fourth 
order Runge Kutta method was compared to an explicit finite- 
difference solution of the partial differential heat equation 
[Eq. (Ill— 36) ] for the same conditions of htating-gas temper- 
ature, propellant properties and variation of h c . Maximum 
difference of 5% between the two solutions for heat-up times 
necessary to raise the surface temperature from 298 to 700°K 
was obtained. 

Equation (III-40) is solved by a standard fourth-order 
Runge-Kutta method simultaneously with the governing equa- 
tions for the chamber flow field to yield the propellant 
surface temperature at any calculated time and position. In 
this way flame spreading rates can be predicted, using the 
ignition criterion as described in Chapter III Section A. 

The non-propellant walls of the port have a coating 
with thermal properties assumed to be equal to those of the 
propellant (see Chapter V, Section A) . Therefore, before 
local ignition T„ s = Tp S . After ignition at a position, 

T ws at that position is calculated continuously by Eq. (HI- 
50). 


The large number of coupled parameters involved in the 
analysis precludes a simple and useful nondimensionalization. 
Moreover, it is doubtful, if any use of normalized parameters 
can provide an easy and straight-forward interpretation of 
the results. Therefore, it was decided to carry on the 
numerical calculations in a dimensional form, using consist- 
ent, metric engineering system of units, thus allowing a 
direct comparison of predictions with experimental measure- 
ments . 
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CHAPTER IV 
NUMERICAL SOLUTION 
A. Implementation Scheme 

Numerical mathematical tec miques were developed for the 
simultaneous integration of the three governing equations, 

Eqs. (III-7), (III-8) , and (III-9) , the two equations for the 
entrance section, Eqs. (III-13) , and (III-14) , and the equa- 
tion, which describes the rate of change of the propellant 
surface temperature, Eq. (III-50) . The numerical steps were 
organized i.ito a program for a large capacity digital computer. 

It can be shown, that the set of governing equations is 
totally hyperbolic in nature . 1" ,190 Indeed, all three eigen- 
values of its characteristic equation are distinct and real. 

The numerical techniques developed aimed at maximum accuracy, 
stable conditions, and computation efficiency. Recently 
obtained experience in solving hyperbolic partial differen- 
tial equations^O was utilized, and a generalized implicit 
scheme, based on central differences in spacewise deriva- 
tives, 191 r 192 was c hosen to solve numerically the governing 
equations. Let the net of points in the t,x-plane be given 
by t = jAt and x = nAx, where j = 0,1,2. . ., and 
n = 0,1,2. . . . The mesh size of the net is determined by 
At and Ax. Then the derivatives of pressure, for instance, 
are expressed in the following difference form: 


iE 

3t 


l£ 

9x 


j+1 j 

P J -p J 
*n H n 

At 


n 


n 


(IV-1) 


6 <*&-*>&> + (1 - 6 ><Pn + l-Pn-l> 


2Ax 


where p^ = p(jAt, nAx). 


The weighting parameter 8 is a real constant, lying 
in the interval 0 £ 6 £ 1. For the implicit formulation no 
stability restriction exist if D.5 £ 9 £ l. 191 Usually a 
value of 0 * 0.6 was used in the calculations. 
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B. Quasi-linearization and Predictor-Corrector Calculations 
of the Non-linear Terms in the Governing Equations 

The governing equations may be written in the following 
matrix form 
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(IV-2) 
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where F [u (t ,x) ,T (t ,x) ,p (t ,x) ] are the functional coefficients 
of the partial spacewise derivatives, and I [u (t ,x) ,T (t ,x) , 
p(t,x)] are the corresponding inhomogeneous terms in Eqs- 
(III-7) , ( III— 8 ) , and (III-9) . 


In order to obtain a system of linear algebraic differ- 
ence equations, which can be solved simultaneously by matrix 
methods, the nonlinear coefficients F(u,T,p) and inhomogen- 
eous terms I(u,T,p) are linearized in a way, described in 
the following paragraphs. 


The inhomogeneous terms were first quasi-linearized, 
according to the scheme :193 
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(IV-3) 


n 


where the partial derivatives, which are actually algebraic 
terms, are evaluated with the flow properties at (j + 0/2, n) , 

.9 . 

such as, for instance, p^ + J ■ p^ + l-(p- J + - p^) . 

n *n 2 '*n r n 


A predictor-corrector technique was applied to the non- 
linear coefficients F(u,T,p), and the partial derivatives in 
Eq. (IV-3) . According to this method, a predictor calcula- 
tion is first made, in which previous-time flow properties 
are used to evaluate all the coefficients F and partial 
derivatives in the quasi-linearized inhomogeneous terms. 

Then, a corrector calculation is made, in which the coeffi- 
cients F and the partial derivatives are evaluated with 
weight-averaged flow properties. The iterations may continue 
on until sufficient convergence is achieved. 


i &*^*W*~ 




* ^ ^ 
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After substitution of the aforementioned finite-differ- 
ence representations into the governing equations, the vel- 
ocity variation equation, Eq. (III-7) , becomes 
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The temperate.® variation equation, Eq. (III-8) , becomes 
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and the pressure variation equations, Eq. (III-9) , becomes 
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The barred parameters are weight-averaged quantities 
as follows: 
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whereas the flow parameters in the partial derivative form 
in Eqs. (IV-4) , (IV-5) , and (IV-6) are 
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The parameters ( ^ 1 )pd are quantities, calculated 
from the latest predictor iteration. In the first step 
predictor calculation they have the values of the previous 
time step. 
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190 194 

As pointed out by Kuo, et al, ' the predictor- 
corrector method, combined with the quasi-linearization 
scheme, is very helpful in stabilizing the numerical solution. 
However, the successful application of this combination to 
the present study does not generalize the approach for all 

similar problems. 1^0 

Negligibly small differences were found between final 
solutions obtained by single-step and multi-step predictor 
calculations. Therefore, to reduce computing time, a single- 
step predictor-corrector calculations are used in the actu- 
al computations. The final solution at every time step is 
obtained implicitly by a matrix method, discussed in a later 
section. 

C. Extraneous Boundary Conditions 

The use of central-difference formulation, as described 
in Eq. (IV-1 ) , for all spacewise derivatives in the governing 
equations, requires six boundary conditions for the solution. 
In other words, the set of finite-difference equations cor- 
responds to a hyperbolic system of 6th order. Therefore, 
three boundary conditions, in addition to the physical ones 
defined by Eqs. (III-13) , (III-14) , and (III-17) , or (III- 
18) are needed. These so-called extraneous boundary condi- 
tions for the hyperbolic system of governing equations stud- 
ied are derived from the compatibility relations at the 
boundaries. The latter are obtained by solving the govern- 
ing equations by the method of characteristics .148 

The compatibility relations along the right-running and 
left-running Mach lines (defined by dx/dt = u ± c, respect- 
ively) , in terms of p, u-characteristics, are 


dt 
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(IV-9 ,IV-10) 


The relation along the particle-path line (defined by 
dx/dt = u) , in terms of T,p-characteristics is 


dT 

dt 


III 


y-1 Tfdpl 
Y p[3tj 


III 


(l - I Zj 

a T v n A 


Y P P 
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In relations (IV-9), (IV-10) , and (IV-11), I u , I?, and 
Ip are the inhomogeneous terms in the governing equations 
[See Eq. (IV-2) ] . 

For a subsonic flow toward the nozzle in the motor which 
is the case during the entire transient in the physical model, 
described in Chapter III, the left-running characteristic line 
is considered at the left boundary, whereas the right -running 
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characteristic line together with the particle-path line are 
considered at the right boundary of the region of numerical 
computation. Taking into account the conditions at the right 
boundary, which is the aft-end of the motor (uniform perime- 
ter, no area change, and no mass addition, see Fig. 2) , the 
compatibility relation along the right-running characteristic, 
Eq. (IV- 9) , becomes 
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The compatibility relation along the left-running 
characteristic, Eq. (IV-10) , is 
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The compatibility relation along the particle-path line, 
Eq. (IV-11) , at the right boundary, becomes 
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Treatment of the Boundary Conditions 


The compatibility relation along the left-running char- 
acteristic line, Eq. (IV-13) , and Eqs. (III-13) and (III-14) 
form a closed system of ordinary differential equations to 
determine the gas velocity, temperature, and pressure at the 
entrance to the propellant section. The system, formed by 
the other two compatibility relations, Eqs. (IV-12) and (IV- 
13) , and Eq. (III-17) , or (III-18) (unchoked, or choked flow, 
respectively) , determines the flow parameters at the chamber 
aft-end. The equations in each system are expressed in a 
finite-difference form and quasi-linearized in an implicit 
form, using predictor-corrector iteration calculation in the 
same method which is utilized for the governing equations. 

As a result, a set of three linear inhomogeneous algebraic 
equations is obtained for every boundary. Each set of equa- 
tions is solved by 3 x 3 matrix, using Kramer's rule, to 
yield the gas velocity, temperature, and pressure at the 
boundary for the next time step calculations. The character- 
istic line segments near the left and right boundary of the 
calculated domain are shown in Figs. 8 and 9, respectively. 
These figures also show the numerical calculation grid in 
an x,t-diagram. 

In Fig. 8 the line sent out from A^ to the left-bound- 
ary point B$,(xg,t3*l) with a slope (dx/dt)jj=u - c is the 
left-running characteristic. The process or determination 
of the boundary values is carried out in the following steps: 


(1) With the calculation completed for the time step . 
j At , all properties are determined at the point 
by linear interpolation between the boundary 
(xi,t3). The position of is taken from the 

previous-time-step calculations, as described 
below. 


(2) The left-boundary system of equations [Eqs. (HI- 

13) , (III-14), and (IV-13)] is solved for u, T, 
and p at B^ (xq , t7 + ^) . The flow parameters, 
present in the coefficie: and inhomogeneous terms 

of these equations, are x. rated with the solutions 
until sufficient convergence is reached. 

(3) The corrected location of A^ at the current time 
is found from the slope of the left-running char- 
acteristic (solid line A^ - Bj, in Fig. 8), i.e., 




(IV-15) 


> 


If l X A2 

repeated with 
subscript "i" 
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0.002, steps (1), ( 2 ), and (3) are 

0.8 x A 2 + 0.2 (x^)^ w ^ ere t ^ ie 

number of the iteration step. 


When the convergence condition for the point A^ is met 
(usually in two iterations) , the distance x of this point 
becomes the distance for the first iteration in the next time 


step, i.e., = x a 2 (dotted vertical line in Fig. 8 ). 

In Fig. 9 the lines sent. out from A^ and A^ to the 
right -boundary point B r (x N ,t3 + l) are the right-running 
characteristic and particle path, respectively. The calcu- 
lation of the right-boundary values by the simultaneous 
solution of the above-mentioned pertinent system of equations 
is similar to that of the left-boundary values. 


The position of the points A^ , A 2 , and A 3 changes very 
slowly with respect to time. The directions of exaggerated 
shifts shown in Figs. 8 and 9 are typical for the chamber- 
filling interval of the transient. 

It should be noted that solutions, based on separate 
integrations along the characteristic lines, 190 have resulted 
in oscillations of the boundary values. The quasi-linearized 
simultaneous-solution method results in a smooth and stable 
solution. The method may be applied, however, with caution 
to other configurations which cause different boundary condi- 
tions, such as zero aft-end velocity during part of the 
transient due to use of a nozzle closure, and change of 
velocity direction at the fore-end due to short igniter 
operation. 

E. Computation Efficiency and Convergence Tests 

The numerical calculation solves Eqs. (IV-4) , (IV-5) , 
and (IV- 6 ), which represent the governing equations in their 
implicit difference algorithm. For the N-l spacewise interi- 
or points considered, 3 (N-l) linear algebraic equations have 
to be solved simultaneously for every time step for 3 (N-l) 
unknowns. The value of N-l depends upon the given lengths 
of the propellant and aft sections, x E - x p , and the spatial 
mesh size. Ax, used in the calculation. The faster the pres- 
surization process, the smaller the value of Ax that must 
be used, and therefore, the larger the number of interior 
calculation points. Also, the time step, or temporal mesh 
size. At, should be kept small enough, so that the variations 
in flow properties with respect to time can be studied, and 
the quasi-linearized nonlinear and inhomogeneous terms can 
be properly represented. 
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For efficient computation, the finite-difference equa- 
tions, Eqs. (IV-4) , (IV-5) , and (IV-6) , are arranged into 

so-called block-tridiagonal matrix form. 190, 195 In this 
particular form, the 3(N-1) x 3(N-1) matrix of the coeffi- 
cients of the unknowns has (3N-5) 3x3 square sub-matrices 

as single elements. An economical solution is obtained 


through block-factorization into two block-bidiagonal 
matrices. 190,195 For accurate computation, the computer 
option of double precision is used in the matrix calculation. 


Variation of the mesh sizes. At and Ax, was carried out 
to test the convergence of the solution. This variation was 
conducted in connection with the so-called Courant-Friedrichs- 
Lewy stability condition for one-step space difference equa- 
tions, which is 

( I u I + c) At/Ax < 1 (IV-16) 


However, when an implicit difference scheme is used, the 
numerical solution becomes unconditionally stable. 191 

It has been found that within the region of tests with 
the experimental model, negligible changes in the solutions 
are obtained when At or Ax are reduced by a factor of 
five. Mesh sizes, for which the maximum value of 
( | u | + c) At/Ax is slightly greater than one, are usually 
used . 


Various tests were performed on the sensitivity of the 
solution to changes in initial conditions and small changes 
in important input parameters, such as the igniter mass flow 
rate, burning rate law, propellant ignition temperature, 
adiabatic flame temperature, and propellant density. In all 
these tests the solutions were bounded and smooth, and changed 
only slightly for small changes in these parameters. This 
demonstrates the existence of neighboring solutions and the 
convergence of the entire numerical solution. 

Appendix D shows a flow chart for the solution of the 
governing equations, and Appendix E shows the overall com- 
puter program flow chart. Lists of the computer program sub- 
routines and the input data cards are given in Appendix F. 

The Fortran IV computer program can be requested from L. H. 
Caveny at the Guggenheim Laboratories of Princeton University. 
Requests for the program should be accompanied by a 9 track 
magnetic tape. Compilation of the program on the IBM 360-91 
requires an area ot 350K for the IBM G level compiler and 
500K for the IBM H level compiler. 
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CHAPTER V 

EXPERIMENTAL INVESTIGATION OF THE STARTING TRANSIENT 


A. The Experimental Moi:or 

An extensive experimental program was carried out to 
provide the empirical data required as input to the theoret- 
ical model and to test the validity of the analysis. A 
laboratory-size solid-propellant window motor was designed 
and manufactured for the experimental investigation. The 
design of the motor was guided by the following primary re- 
quirements. It is necessary that: 

(1) The experimental set-up be as close as possible to 
a practical rocket motor, both in configuration 
and range of operation; 

(2) The motor has sufficient flexibility to permit 
wide-range diagnostic studies involving measure- 
ments of different types; 

(3) The analytical model and the experimental model be 
compatible for tests of the theory without unnec- 
essary complications, and 

(4) The size of the motor be easily manageable and 
consistent with laboratory testing procedures. 

Schematic longitudinal and cross sections of the exper- 
imental motor are shown in Figs. 2 and 3, respectively. 

Figure 10 shows an assembly drawing and Fig. 11 is an explod- 
ed view of the laboratory-size rectangular window motor. 

The propellant slabs were cast into brass trays. To insure 
good propellant-to-tray bonds, the inner surface of the tray 
was sand-blasted and coated with a thin layer of Duro-Plastic 
Epoxe glue (available from Woodhill Chemical Sales Corp.). 

The exposed propellant surface was cast and polymerized a- 
gainst the polished surface of a Teflon strip. Several types 
of trays with two propellant web thicknesses (0.635 and 0.318 
cm) were used. Each tray had a constant web thickness. The 
two trays used at a test are so positioned in the propellant 
section, as to create an initial uniform port of 2.54 by 
0.635 cm along the entire motor, from the end of the igniter 
nozzle to the entrance to the motor nozzle. The primary 
feature of the slab-grain geometry is that it permits direct 
photographic observations during ignition and flame spreading, 
and at the same time the chamber conditions realistically 
correspond to the conditions in a practical rocket motor. 
Furthermore, many diagnostic studies, such as pressure, tem- 
perature and heat-flux measurements can be conveniently con- 
ducted . 
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One side wall of the port is a two-part Plexiglas window. 
Tests with fast water quench revealed that no Plexiglas abla- 
tion occurs prior to the time of peak pressure up to 100 msec. 
The inner part of the window was replaced for each test. The 
window is sealed by an o-ring gasket and retained by a cover 
connected to the motor block by 26 high-strength bolts. For 
tests in which measurements other than photographic were made, 
the window was replaced by a stainless steel part having the 
same dimensions. 

The side wall of the port channel opposite the window 
is used for the instrumentation ports. Five Dynisco PT76 high 
frequency pressure transducers can be installed along the port 
12.7 cm apart, starting from the propellant leading edge. A 
burst diaphragm is connected to each end of the motor for over 
pressure protection. Three additional ports are available for 
thermocouple probes and heat-flux gauges. The total ten ports 
can be used in different combinations for techniques, such as 
depressurization, liquid quenching, side ignition, and purging 

The motor nozzle, made of copper, consists of a rectangu- 
lar convergent section and a precisely defined rectangular 
throat. Four nozzles were manufactured for series of tests 
with port-to-throat area ratio of 1.06, 1.2, 1.5 and 2.0. 

There was no measurable erosion in all tests conducted. The 
nozzle block is bolted to the aft end of the motor with a 
connector, which leads the combustion gases from the nozzle 
to the exhaust vent system. 

To reduce the heat loss to the chamber walls and, thus, 
to make the experimental motor more compatible with operation- 
al motors, the entrance and aft sections of the motor and the 
nonpropellant walls of the port (except the Plexiglas window, 
when used) were coated before each test with a thin polymeric 
layer consisting of a mixture of 50% PBAA/EPON binder and 50% 
TiC>2 powder. This layer ablated only slightly during the 
starting transient. In the analytical model it is assumed 
that this layer has the same thermal properties as the solid 
propellant. A special test (described in Section B, Chapter 
VI) was conducted to verify the validity of this assumption. 

B. The Ignition System 

The igniter source used in this study was a head -end 
pyrogen igniter, which burned gaseous methane and oxygen in 
a controllable manner. The system is capable of burning a 
wide variety of gaseous oxidizers and hydrocarbon fuels (or 
mixtures of gases) and thus operating with a wide range of 
igniter temperatures. The methane and oxygen were of research 
purity and supplied by the Matheso- Gas Products Co. 

The specially-designed igniter combustion chamber achiev- 
ed 80% of the equilibrium operating pressure within 3 msec. 

The chamber consists of a small stainless steel tube (0.9 cm 
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in diameter and 7 cm .long) , placed in a copper block. The 
injector head comprises two carefully calibrated unlike 
impingement doublets. Fast ignition is achieved by four 
small 10 mm spark plugs (type Z-6, produced by Champion 
Spark Plug Co.), positioned along the igniter chamber. A 
small burst disk is also installed in the system for over- 
pressure protection. The pressure in the combustion chamber 
was measured by a Dynisco PT76 pressure transducer. The ig- 
niter gases were discharged into the entrance section of the 
motor through a convergent nozzle. 

A controllable feed system supplies the gaseous reac- 
tants to the igniter combustion chamber. Each reactant, 
coming from the high-pressure cylinder, passes through a 
dome pressure regulator (Model 15H, 0-2000 psi, made by Grove 
Reg. Co.), which has remotely preset outlet pressure. At 
this pressure, the gas passes through a check valve (Type 
K-1359-8-S made by Kohler Co.) and an on-off pneumatically- 
operated valve (Model PV20F, made by Marotta Valve Corp.) 
before entering the injector. The two on-off valves on the 
fend lines are operated simultaneously by a single solenoid 
valve (Marotta, Model PV74) . The total opening time of the 
on-off valves is 2 msec after the stem moves off the seat. 

The capacitance of the feed lines between the on-off valves 
and the injector was minimized to assure fast start-up and 
shut-off of the igniter. The feed pressures just upstream 
of the injector were measured by Dynisco PT76 pressure trans- 
ducers. A cold-flow adjustment of the feed pressures was 
performed before each run. 

A purging subsystem is an integral part of the control- 
lable igniter. It consists of a nitrogen supply line, a 
solenoid valve (Marotta, model PV74) and two check valves 
(Kohler, type K-1359-2) . The nitrogen purge enters the 
system just downstream of the on-off valves. The igniter 
chamber and rocket motor were flushed with nitrogen before 
and following each test. 

Figure 12 depicts the assembled igniter combustion 
chamber, injector, two on-off valves, fuel-line pressure 
transducer, fuel-line purge check valve, and the sonic igni- 
ter nozzle. The entire assembled ignition system on the test 
bench is shown in Fig. 13. 

For all runs the fuel and oxidizer feed pressures were 
selected to be higher than the maximum pressure in the motor. 
In this way back flow from the motor to the ignition system 
and erosion of the injector orifices were prevented. Com- 
binations of injector orifices and igniter nozzles were used 
to allow experiments with different igniter mass flow rates 
and equivalence ratios. The equilibrium pressure in the ig- 
niter combustion chamber during the induction period was 
usually between 20 and 30 atm. The injector orifices 
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were always choked during the induction interval and served 
as metering orifices. Real-gas effects f important at the 
high feeding pressures, were accounted for in the precise 
determination of the igniter mass flow rates. 1^6 The 
discharge coefficients of the various injection orifices 
were between 0.81 for the smallest fuel orifices and 0.92 
for the largest oxidizer orifices. The c* efficiency of 
the igniter chamber (defined as the ratio of c* calculated 
from the measured mass flow rate and pressure in the combus- 
tion chamber and c* calculated from adiabatic thermochemis- 
try) was 95%. 

As part of the efforts to simulate real rocket motors 
having pyrogen igniters, the ignition system operated during 
the entire starting transient with constant feed pressures. 
During the last part of the flame spreading period and there- 
after, the igniter nozzle becomes unchoked, and the pressure 
in the igniter combustion chamber increases in accordance 
with the rising pressure in the motor. The quasi-steady 
igniter mass flow rate after the unchokina ' s somewhat less 
than the previous steady-state value, acr ding to the par- 
ticular combination of pressures. Sincf :he ef^ -ct of the 
igniter mass flow rate during the last -wages of the tran- 
sient is very small (as compared with that of mass addition 
from the burning propellant) , it is assumed in the theoretical 
model that during this period the igniter operates with the 
same mass flow rate as before the unchoking. 

Figure 14 shows a typical schematic variation with time 
of: (a) the quasi-steady igniter mass flow rate calculated 

from experimental pressure-time measurements, (b) the igni- 
ter mass flow rate used in the numerical calculations, and 
(c) the measured head-end motor pressure. 

The sequence of events in the operation of the ignition 
system starts with establishing the desired feed pressures 
by adjusting remotely the two dome pressure regulators. An 
automatic electrical sequential cam system actuates the 
following operations. When the firing circuit is energized, 
the automatic sequencer first activates the spark plugs. 

Then the solenoid valve (which controls the feed valves) 
opens and the on-off feed valves unclose to start the flow 
of the reactants into the igniter combustion chamber. The 
igniter action is terminated when the sequencer signals the 
valves to close and afterwards power to the spark plugs is 
turned off. The duration of the igniter action can be alter- 
ed simply by adjusting the timing of the shut-off cam. The 
easy variation of duration time, igniter-gas temperature and 
composition, and igniter mass flow rate are the main advan- 
tages of the gaseous ignition system. 

The gaseous pyrogen igniter was usually operated slight- 
ly on the fuel-rich side (with an equivalence ratio ranging 
between 1.02 and 1.08). The stoichiometric mixture ratio of 
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the oxygen-methane system is 4:1. Thermochemical calcula- j 

tions^o yield at 22 atm an adiabatic flame temperature of j 

3460°K and major product mole fractions of 43% H2O, 15% CO ; 

and 13% C02* f 

I 

For the same fuel-oxidizer combination, the adiabatic 
flame temperature of the gaseous igniter does not change 
significantly when changing the equivalence ratio between 
1.0 and 1.4. Therefore, wide range of igniter gas tempera- 
tures may be obtained in experimental tests by using differ- 
ent fuel-oxidizer combinations. 

Figures 15 and 16 show the experimental igniter-motor 
system mounted on the test bench in the Solid Propellant 
Test Cell. 

C. Instrumentation and Data Reduction 

Two types of physical measurements were conducted in 
the course of the experimental research program: pressure 

measurements and photographic observations. In each test, 
the fuel and oxidizer feed pressures (50 to 80 atm) and the 
pressure in the igniter combustion chamber were measured. 

In addition, up to five pressure measurements along the cham- 
ber port we re made. The Dymsco pressure transducers used 
are of the bonded strain gauge type. Model PT76, suitable 
for dynamic pressure measurements. The natural frequency of 
these transducers as stated by the manufacturer is 25,000 
cps. From the investigations of Thomas and Laytonl*9-201 it 
was deduced that the amplitude frequency response is flat up 
to 8,000 cps. The frequency range of the transient pressure 
phenomenon investigated in this study is 40 to 400 cps. The 
diaphragm of the transducers used lor pressure measurements 
in the igniter chamber and in the motor was protected from 
the hot combustion gases by a 1/8 inch thick disk made of 
General Electric RTV-580, which was bonded to the diaphragm 
by RTV-108 adhesive. Shock-tube experiments showed that the 
natural frequency of the protected PT - ^ transducer is 24,000 
cps. The installation of the transducers was performed ac- 
cording to the research experience accumulated in the AMS 
department of Princeton Universityl*9“203 and elsewhere. 

205 Figure 17 depicts the back side of the experimental 
motor showing the f've pressure transducers installed along 
the motor. 

The transducers in use were calibrated in position before 
each experimental test. The calibration was performed by 
pressurizing the assembled rocket motor with nitrogen gas and 
by regulating the nitrogen pressure. A master pressure gauge \ 

(Heise, Type H8755R, 0 to 2000 psi) was used for the calibra- \ 

tion in steps of 100 psi. The master gauge was periodically 
checked and calibrated with an Amthor (Type 472, 0 to 2000 
psi) dead-weight tester. Also periodically a calibration of 
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the pressure transducers after a test was performed with the 
transducers still in position. Figure 38 shows the Solid 
Propellant Test Cell Calibration and Firing Console and part 
of the Instrumentation Console. 

The /pressure transducers were driven with B&F Input 
Signal Conditioners. These units have a ripple level of less 
than 0.5' milli'* jlts peak to peak and were used to supply 10 
volts d.c. excitation signal. The pressure signals coming 
from the transducers were amplified by DANA (Models No. 3420 
and 3520) d.c. amplifiers. Model 3520 has built-in filters 
with band pass widths ranging from wide band and 10 KC down 
to 0.01 KC. The amplified pressure signals were recorded on 
a Honeywell Visicorder (Model No. 1508) dir jct-recording 
oscillograph. The Visicorder used is capable of recording 
simultaneously on up to 12 active channels, using light- 
sensitive paper. The galvanometers used were Honeywell 
Series "M" sub-miniature, fluid damped, Type No. M1000 and 
M1650 with nominal undamped natural frequency of 1000 and 
1650 cps, respectively. The visicorder was usually operated 
at its maximum speed, 120 inches of paper per second, with 
10 millisecond timing marks. More details on the instrumen- 
tation used for the pressure measurements and recording can 
be found in the Instrumentation Manual206 f or the Solid- 
Propellant Test Cells in the Guggenheim Laboratories of 
Princeton University, and in Ref. 74. 


The photographic observations were made by taking high- 
speed motion pictures with a 16 mm Wollensak Fastax camera 
(Model WF-16) . Both black and white (using DuPont 931-A 
rapid reversal film) and color (using Kodak Ektachrome EF 
film) movies were taken. Operating speeds of 1200, 2000 
and 4000 fps were used. The camera was synchronized with 
the experimental test so that the operating speed of the 
camera was attained before first ignition. Two neon timing 
lamps were used with the photography. The synchronization 
timing lamp was connected to a cam of the ignition sequence 
timer to generate a signal shortly after ignition, which was 
recorded simultaneously on the oscillograph and the film for 
synchronization. The second lamp was connected to a signal 
generator and produces timing m«rks every millisecond of real 
time on the film. 


Data reduction of eight simultaneous records per test 
is a tremendous ->b if not performed in an automatic manner. 

In the later stages of the study analog-to-digital conversion 
(A - D) was used to reduce and process the data. The Daiacom 
^Data Acquisition System (Model 8015) used in this study is 
capable of acquiring analog data from up to 64 sources at a 
^tal sampling rate of 40,000 per sec. The analog channels 
multiplexed, converted to 12-bit. binary words by a Ray- 
thefe^n A-D Converter and recorded on a 9 -track magnetic tape 
compatible with the IBM 360/91 computer. The data is recorded 
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in blocks which are stored in the system memory with a capac- 
ity of 2048 bytes. When a block is complete it is transfer- 
red automatically to the tape, while the next block is being 
filled. The Datacom system was used simultaneously with the 
Visicorder oscillograph during both calibration and test re- 
cording. Proper wiring and grounding was made to che Datacom 
system to make it compatible with the Visicorder. The data 
recorded on the tape was demultiplexed, converted back to 
pressure readings, and analyzed using an IBM 360/91 computer 
by means of a program specially written for this purpose. 
Typical use of the Datacom system in this study included re- 
cording of the five motor-pressure measurements at total 
sampling rate of 10,000 readings per second (or 2,000 readings 
per channel per sec) and block length of 1800 bytes with the 
pressurization part of the transient included in one block. 

A mean deviation of less than 1% was obtained between pressure 
vs time plots from manual reduction of Visicorder records and 
those processed by using the Datacom system. The mean devi- 
ation is defined here as the ratio between the area enclosed 
by the twc p - t plots obtained from the two systems and 
the area enclosed by the Visicorder plot and the time axis. 

Figure 19 shows the Solid Propellant Test Cell Instru- 
mentation Console and the Datacom Data Acquisition System 
with the related wiring before a test. 

D. Experimental Test Series 

The major parameters varied in the experimental study 
were the port-to-throat area ratio and the igniter mass flow 
rate. According to the four port-to-throat aree ratios used, 
2.0, 1.5, 1.2 and 1.06, four series of runs were conducted, 
designated C, D, E and F, respectively. In each series 5 to 
6 tests were carried out with different igniter mass flow 
rates, ranging from 8.5 to 21.5 g/sec. In all these runs the 
gaseous igniter was operated at close to stoichiometric fuel- 
oxidizer mixture ratio, yielding a theoretical adiabatic 
flame temperature of 3460°K at 21 atm as described in Section 
B of this Chapter. The same composite propellant, the prop- 
erties of which are listed in Table 1, was used in all experi- 
mental tests. The purpose of these run series was to test 
the entire experimental set-up, investigate the effects of 
port-to-throat area ratio and igniter mass flow rate on the 
various characteristic transient parameters of interest and 
check the validity of the analytical model. Several early 
experiments were conducted at the start of the experimental 
program to check reproducibility of test results. For the 
same igniter mass flow rate very good reproducibility of 
pressurization rates and maximum pressure was obtained. In 
all test motor runs p(t,x) measurements were taken, and in 
some runs photographic observations were made. In addition 
to the aforementioned four series of runs several particular- 
ly diagnostic experimental tests were carried out as described 
in the following paragraphs. 
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An experimental test was conducted at Ap/A t = 1.2 with 
the igniter operating at an equivalence ratio of 1.35, which 
yields a theoretical adiabatic flame temperature of 3360°K, 
to confirm experimentally the findings of the theoretical 
calculations about the effect of igniter gas temperature. 

An experimental firing was carried out with shorter pro- 
pellant slabs (33.0 cm long) in the upstream part of the pro- 
pellant section. The PBAA/Ti02 mixture v/as applied in a thin 
layer over the downstream tray surface to check: (1) the 

validity of the assumption that PBAA/TiC >2 layer has approx- 
imately the same thermal properties as the propellant, and 
(2) the effect of propellant grain length for the same config- 
uration on the transient. 

Liquid quench experiments were conducted to find the 
erosive burning exponent of the propellant and to check the 
validity of the analytical model in the post-transient period 
by measurements of distance burned. A quench system based on 
the study of Strand and Gerber^O? was designed, built, and 
successfully operated. Water was injected simultaneously 
through five ports along the motor. Check valves just up- 
stream of the injector heads prevented entry of the combus- 
tion gases into the system. The water jets bounced oi- the 
opposite wall in small droplets and thus propellant surface 
cutting was avoided. Small amounts of water (as little as 
5 gms) injected in 10 msec caused a complete extinguishment 
in 3 msec. It was discovered in the quench experiments that 
the Plexiglas window does not ablate even when the test time 
is extended to the period of steady-state operation. This 
means that no corrections are necessary in the analytical 
model for port area changes and for mass addition resulting 
from Plexiglas regression during the transient period. 

Techniques for direct photographic measurements of burn- 
ing rate during the transient and thereafter were attempted 
during the course of this study. The techniques were not 
developed fully since burning rate measurements were not con- 
sidered a primary goal in the study. A promising technique 
tried with a partial success in the investigation will be 
mentioned. The propellant was cast simultaneously in the tv/o 
opposing trays and then bonded to and poiymerizt J wit .i a 
Plexiglas window, which was connected by screws to the trays. 
The propellant grain-window assembly is depicted ir. Fig. 20. 
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CHAPTER VI 

DISCUSSION OF RESULTS 


Diagnostic experimental tests and corresponding numeri- 
cal calculations weie conducted with AP composite propellant 
at port-to-throat area ratios of 2.0, 1.5, 1.2 and 1.06 with 
igniter mass flow rates ranging from 8.5 to 20.5 g/sec. 

Table 1 lists the baseline values of propellant properties 
and motor parameters used in the calculations. These values 
correspond closely to the experimental condi .ions in all 
tests. A datum case was defined at Ap/A^- = 1.2 and mig = 

15.0 g/sec using the baseline values of Table 1. It was so 
selected because port-to-throat area ratio of 1.2 and oper- 
ating pressure of about 40 atm represent a realistic goal 
for high-performance motors that designers presently cannot 
analyze with confidence. The results of the experimental 
study are listed in Table 2. Table 3 summaj.izes the results 
of the theoretical parametric study about the datum case and 
calculations corresponding to experimental tests. A list of 
input values for the datum case calculations is presented 
in Table Fi (see end cf Appendix F) . Typical calculated 
Mach and Reynolds numbers during the transient for the dif- 
ferent port-to-throat area ratios tested are given in Table 4. 

A. Comparison Between the Theoretically Predicted and 
Experimentally Measured Pressure-Time-Space Traces 

Typical starting transient p(t,x) measurements, ob- 
tained during expe: imental tests at the five measurement 
stations (denoted in Fig. 2) at port-to-throat area ratios 
c" 2.0, 1.5, 1.2, and 1.06 are sh^v.’n ir. Figs. 21 through 24, 
rc pectively. By recalling Fig. 1, it is seen that HVT motors 
c e operating in a pressurization regime which is totally 
different from that in the ignition studies of Most and Sum- 
merfield.5 The plots show clearly the strong time and space- 
dependence of motor pressure during the transient and thus 
demonstrate the necessity of this analysis. Consider the 
starting pressure transient depicted in Fig. 23. Following 
a small pressure wave caused by the onset of igniter flow, 
ail almost uniform pressure distribution at 1.2 atm was estab- 
lished in the motor cavity. At 36 msec after the igniter 
start, a flame was observed (on high speed photographs) at 
the leading edges of the propexlant slabs and flame spread- 
ing begins, accompanied by increase of pressure. The pressure- 
rise at station 3 (midd.1e of motor) starts later but increases 
at the fastest rate. Pressure peaks are achieved at 72 msec 
at the fore-end (station 1) and at 73.5 msec at the aft-end 
(station 5) . 

Axial variation of pressure increases as Ap/A t decreases. 
In all tests, the pressure measured at station 5 was almost 
identical to that measured at station 1 (fore-end of motor) 
due to the igniter-motor configuration and the low Mach num- 
bers at the fore-end part of the propellant section. 
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Comparisons between typical measured and calculated (by 
the analytical model) pressure- time curves at stations 1, 3, 
and 5 for port-to-throat area ratios of 1.5, 1.2, and 1.06 
are shown in Pigs. 25, 26, and 27, respectively. Also shown 
in these figures are the pressure-time traces calculated for 
no erosive burning and the calculated pressure-time point of 
flame spreading completion. In general, very good agreement 
between measured and calculated pressure-time plots is observ- 
ed. The only input variable adjusted to achieve agreement 
was the erosive burning exponent 6. However, in future 6 
can be established a priori by erosive burning tests (such 
as burnt distance vs x measurements after extinguishment 
shortly after the transient) . It is emphasized that no ad- 
justments in parameters are made during the parametric stud- 
ies. The solution of the analytical model predicts precisely 
the initial development of pressure distribution at the onset 
of igniter and the following nearly constant pressure distri- 
bution during the induction interval (not shown in Figs. 25 
through 27) . Agreement within 10% between measured and cal- 
culated time to first ignition (i.e., occurrence of first 
flame on propellant surface) and time to attain the maximum 
fore-end pressure is obtained for all comparisons at port- 
to-throat area ratio of 1.2, 1.5, and 2.0. No ignition time 
measurements were made at Ap/A^ = 1.06. The agreement between 
measured and calculated time to attain the maximum fore-end 
pressure for this case is within 15%. 

An agreement within 10% (in most cases less than 5%) is 
obtained between calculated and measured pressure-peak values 
at all stations and for all port-to-throat area ratios inves- 
tigated. The largest disagreement is observed at A p /A t = 2.0. 

For A p /A t = 1.2, 1.5 and 2.0, during the flame spreading 
period the p calculated pressure is usually slightly lower than 
the measured pressure, whereas during the chamber filling, 
the calculated pressure is slightly higher than the measured 
pressure. It should be mentioned that this type of deviation 
between the calculated and measured pressure transient is the 
same as the one obtained by Most and Summerfield^ *15 for 
motors with high port-to-throat area ratio using the dynamic- 
temperature homogeneous model. A possible explanation for this 
deviation may be, for instance, increased heat transfer to the 
propellant due to surface roughness at the early stages of 
flame spreading, and attenuated heat transfer tc the unignited 
propellant surface at the later stages of flame spreading due 
to considerable blowing from the surface, when the ignition 
temperature is approached. 

For Ap/A c - 1.06, all calculated pressure-time plots pre- 
cede the experimental ones. The maximum deviation in the 
time-axis direction for the range of igniter mass flow rates 
tested is 15%. For the same igniter mass flow rate, as the 
port-to-throat area decreases the velocity over the propellant 
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leading edge increases and the pressure decreases. The phe- 
nomenon of delayed ignition at high gas velocities was observ- 
ed and discussed by McCune,29 Kling, et al, ^ Keller, Baer, 
and Ryan, 6 3, 64 anc j Bastress and Niessen.^ 9 The explanation 
of Keller, Baer, and Ryan is recalled here; as the gas veloc- 
ity is increased, reactive species at or near the surface are 
being increasingly diluted and swept away by the fast-moving 
gas. This is consistent with the research of Kashiwagi and 
co-workers , 60 which shows that the gas-phase kinetic processes 
are important during ignition, and increased flow velocity 
shifts the ignition position downstream in the case of a gas- 
phase ignition model. At Ap/A^ = 1.06 and low to moderate 
igniter mass flow rates the pressure in the experimental motor 
(cs i - every HVT motor) is very close to the ambient The 
increase of ignition delay time when lowering the pressure 
(and approaching the atmospheric pressure; at moderate to high 
surface heat fluxes has been observed and studied by many in- 
vestigators. 32 , 34 , 69 , 70 ,208 , 209 j n this study, typical cal- 
culated convective heat flux at the propellant leading edge 
(average for the induction interval) for Ap/A t = 1.06 and 
mig = 17.2 g/sec was 42 cal/cm2-sec. 

It should be noted that in many tests using photography 
start of pressurization in the motor chamber was measured 
before any ignition has been observed (typically a pressure- 
rise of 1 to 2 atm over the pre-ignition level) . This may 
be a result of a pre-ignition decomposition-gasification of 
the solid-propellant. The reduced cooling of the igniter 
gas due to surface temperature increase, as calculated by 
the analytical model, does not cause a substantial change in 
chamber pressure. 

The contribution of erosive burning to the increase of 
pressure is seen in Figs. 25 through 27. It may be defined 
as the ratio of maximum fore-end pressure calculated with 
erosive burning to the maximum fore-end pressure calculated 
without erosive burning. For the propellant used in this 
study (see Table 1) and the burning rate law given by Eq. 
(III-35) , the aforementioned ratio is approximately the same 
for all values of A p /At tested, and equals 1.5. This means 
that in terms of maximum pressure, the burning rate increase 
along the port due to erosive burning is fully compensated 
by the burning rate decrease due to the significant pressure 
drop. Indeed, the high value of the pressure exponent (n = 
0.4) shows high pressure sensitivity, whereas the large value 
of the erosive burning exponent (8 = 105) represents low 
sensitivity to erosive burning. With respect to motor design, 
this conclusion points out the feasibility of using propel- 
lants with relatively large values of n and B in high- 
performance rocket motors without facing unacceptably high 
fore-end pressures. In the range tested, the erosive burning 
contribution is only slightly affected by the igniter mass 
flow rate, due to its small magnitude as compared to the mass 
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added from the propellant burning. The ratio of steady-state 
igniter mass flow rate (m-j a ) to maximum mass flow rate out of 
the nozzle varied between 3 and 6%, as listed in Table 3. 

One of the principal advantages of the analytical model 
is the ability to calculate spatial, as well as temporal, 
variations of gas parameters. Measured and calculated spatial 
pressure distributions at four illustrative times (correspond- 
ing to different stages of the transient) of Run D-9 (A p /A t = 
1.5) are compared in Fig. 28. Calculated spatidl velocity 
and temperature distributions at the same times arc plotted 
in Figs. 29 and 30, respectively. These plots correspond to 
the pressure-time traces in Figs. 22 and 25. The changes in 
velocity distribution during the transient are well depict ;d 
in Fig. 29. During the induction interval there is a 1~. je 
velocity decrease along the port due to gas cooling. Du. ing 
flame spreading there is velocity increase upstream of the 
ignition front (due to mass addition) and velocity decrease 
downstream of the front. After the completion of flame 
spreading the gas velocity increases considerably in the 
propellant section and changes only slightly in the aft sec- 
tion. 


To show the effect of port-to-throat area ratio, the 
pressure, velocity, and temperature distributions at the 
different stages of the transient for a test with Ap/A^ - 
1.06 are plotted in Figs. 31, 32, and 33, respectively T 
These plots correspond to the pressure-time traces in Figs. 

24 and 27. 

The significant changes in gas parameter distributions 
during the transient are well described in Figs. 28 through 
33 and they emphasize the need for a distance and time 
(p(t,x)] analysis for proper prediction and control of the 
starting transient of HVT motors. 

B. Ignition Delay Times, Flame Spreading, and Pressurization 
Rates 


Most of the photographic observations were made at port- 
to-throat area ratio of 1.2. The event of first ignition was 
relatively easy to detect, whereas flame spreading photography 
was only partially successful as explained later in this sec- 
tion. First ignition, when detected, always appeared at the 
propellant leading edge. Figure 34 shows the measured and 
calculated effect of igniter mass flow rate on the ignition 
delay time t^g (defined as the time from onset of igniter 
to the first appearance of flame on the propellant surface) 
for motor test series "E" (Ap/A t = 1.2). From the thermal 
ignition theory, assuming an F inert propellant being heated 
by a constant heat flux at the surface and a prescribed sur- 
face temperature for ignition, it is deduced that the square 
root of the ignition delay time is inversely proportional to 


-69- 


; \ 


the surface heat flux. If a turbulent flew correlation for 
the convective heat-transfer coefficient Is used [such as 
Dittus-Boelter correlation, Eq. (III-21)], it follows that 
tfg ^ m-Lg -1 * 6 . However, actually, the surface heat flux is 
not constant during the induction interval, because m^- is 
not a step function of time (see Fig. 14) , the heat- transfer 
coefficient is coupled to the varying gas dynamics in the 
motor chamber, and the rising surface temperature decreases 
the driving temperature difference. Figure 35 shows the 
variation with time of the calculated heat flux at the pro- 
pellant surface at three different positions along the motor 
corresponding to pressure measuring stations 1, 3, and 5 (see 
Fig. 2) for a test with Ap/A^. = 1.2. The calculation is made 
in accordance with the experimental variation of igniter mass 
flow rate (see Fig. 14) , such that steady-state value and 80% 
of it are obtained 12 msec and 2 msec after the onset of ig- 
niter, respectively. The shape of heat-flux variations shown 
in Fig. 35 are typical for the calculated distributions in 
this study. Figure 36 shows the calculated propellant sur- 
face temperature vs time at the same positions and for the 
heat flux variations as shown in Fig. 35. 

The slope of the log t^g vs log mig plot on Fig. 34 is 
-1.5. Similar plot of three measured ignition times (not 
shown here) for Ap/A t - 2.0 (see Table 2) shows a slope of 
-1.42. For constant q and m^g, it would mean that tur- 
bulent flow (although not fully developed) existed over the 
propellant leading edge. Calculated ignition delay times vs 
igniter mass flow rate for the four port-to-throat area 
ratios tested are shown in Fig. 37. The slope of the calcu- 
lated lines varies between -1.25 and -1.35, which means that 
the aforementioned factors (accounted for in the analysis) 
affect the calculated ignition delay times. For the same 
igniter mass flow rate in the range investigated, calculated 
ignition delay times decrease slightly when increasing A p /A t 
between 1.06 and 1.5. No photographic measurements were^ 
made at A p /A t = 1.06. 

High-speed photography of fast flame spreading is a 
difficult task. Flame propagation appears to be more like 
an increasing cross-sectional density of ignited points rather 
than a well-defined advancing flame front. Ignition spots 
appear at random downstream of the fully ignited surface. 

This may be caused by an increased heat transfer at local 
surface irregularities. The same phenomenon was observed by 
Mitchell and Ryan, 31 ' 8 ^ by Keller, Baer, and Ryan, 64 by 
McAlevy, et al , ® and by Most and Summer fie Id .*» 1* Figure 
38 shows calculated and measured locations of the flame 
front versus time for an experimental test with port-to- 
throat area ratio of 1.2. Typically, at the beginning of 
flame propagation experimentally measured rates are higher 
than the calculated and this difference reverses toward the 
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920 cm/ sec to 19,800 cm/sec at the start and the end of the 
process, respectively, were obtained for Ap/A t = 2.0 and 
mj_g = 19.5 g/sec. Calculated initial flame spreading rates 
increase substantially with increasing igniter mass flow rate 
but do not show an effect of Ap/At, as seen in Fig. 39. The 
few experimental data are in agreement with the former obser- 
vation but are insufficient to verify the latter. Calculated 
flame spreading times, however, show that for the same igniter 
mass flow rate (and chamber geometry) , the average flame 
spreading rate decreases with Ap/A t . The larger throat area 
results in lower pressures (and thus lower burning rates) and 
lower heat flux, because for the propellant and geometry used 
in this study, the effect of lower pressure overcomes the 
effect of higher velocities. 

With respect to the total pressure rise, the calculated 
pressures at the instant of flame spreading completion vary 
from 25 to 30% of the maximum pressure for A p /At = 2.0 to 
50 to 55% for Ap/A t = 1.06. 

For all port-to-throat area ratios tested, the highest 
pressurization rates in the tests were measured at station 3 
(middle of motor), as seen in Table 2. The maximum pressur- 
ization rates at the motor fore-end (station 1) were slightly 
lower. As far as the numerical calculations are concerned, 
the maximum calculated pressurization rates were usually ob- 
tained at a position between stations 2 and 3. Both measured 
and calculated maximum pressurization rates are obtained at 
the instant of flame spreading completion Dr immediately 
thereafter. 

Measured and calculated maximum pressurization rates at 
station 1 vs igniter mass flow rate for the different Ap/A t 
rested are plotted in Figs. 40 and 41, respectively. Fore- 
“*nd maximum pressurization rates as low as 3,000 atm/ sec for 
Ap/Af * 1-06 and m^g = 15.1 g/sec, and as high as 6,800 atm/ 
sec for Ag/A-t- = 2.0 and m^g = 19.5 g/sec were measured. The 
calculated rates are greater than the measured ones on the 
average by a factor of 2. The maximum pressurization rate 
increases with both the igniter mass flow rate and the port- 
to-throat area ratio. The former effect is less pronounced 
experimentally than theoretically. The maximum pressurization 
rates decrease as the distance from the fore-end increases 
and this decrease is larger as Ap/A*. is smaller. 

The test with shorter propellant trays (see Table 2) 
demonstrates the validity of the assumption that the PBAA/ 

TiC >2 coating and Plexiglas window can be treated as having 
the same thermal properties as the propellant. It also shows 
the effect of burning-to-throat area ratio on pressurization 
rates and maximum pressures. 



C. The Thrust Transient 


Prediction and control of the thrust transient require 
knowledge of the transient stagnation pressure at the nozzle 
entrance, Pas' besides qualitative knowledge of the unsteady 
phenomena taking place in the nozzle. The calculated tran- 
sient variation with time of the fore-end and aft-end (nozzle 
entrance) static pressures and the aft-end stagnation pressure 
for a test with A D /A t = 1.2 are shown in the upper part of 
Fig. 42. The lower part of the same figure shows the varia- 
tion with time during the transient of the ratio of aft-end 
stagnation pressure (p a s) to fore-end static pressure (p]_) . 

The stagnation pressure is calculated from the static pres- 
sure and velocity. During the induction interval p as is 
higher than p^ by 15 to 20% since the effect of gas cooling 
along the port is greater than the effect of friction. After 
first ignition, the ratio Pas/Pl drops gradually to a mini- 
mum of 0.73 and increases afterward to the steady-state value 
of 0.844. After the completion of flame spreading the stag- 
nation pressure decreases along the port due to the dominant 
effect of mass addition. Figure 42 illustrates the conclu- 
sion that calculation of the thrust transient of HVT motors 
by using measured or calculated fore-end static (or any 
"uniform") pressure is improper and misleading. 21 ® 

Figure 43 shows the effect of Ap/A t on calculated (in 
this study) ratios of maximum static and stagnation pressures 
at motor aft-end to maximum motor fore-end pressure. The 
obtained values correspond closely to the values calculated 
by using the steady-state approach of Price 21 ® and the base- 
line quantities listed in Table 1. 

D. Parametric Studies 


1. Effect of Igniter Mass Flow Rate 

The igniter mass flow rate was the major variable 
in the experimental study. Its large effect on the ignition 
delay time and flame spreading rates, and small effect on the 
maximum pressurization rates is described in Section B of 
this Chapter. Experimental and calculated maximum fore-end 
(station 1) pressures vs igniter mass flow rate for the dif- 
ferent port-to-throat area ratios tested are plotted in Fig. 
44. A cross-plot of calculated maximum fore-end pressures 
vs Ap/At for two igniter mass flow rates is shown in Fig. 45. 
The maximum fore-end (and in other stations as seen from 
Tables 2 and 3) pressure increases slightly when increasing 
mig, as expected. The calculated rate of increase is higher 
at Ap/At = 2.0, than at the other port-to-throat area ratios. 
For a given igniter mass flow rate (and the constant config- 
uration under consideration) , the pressure level rises fast 
when increasing Ap/At (by decreasing the nozzle throat area) , 
as illustrated in Fig. 45. 
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The aforementioned effects and trends apply to the range 
of m^g tested. They may change direction or be overcome by 
other effects at extreme cases of marginal m^ (causing a 
hangfire) or very high m^g. An experimental test at Ap/A*- 
= 1.2 and m^g = 8.5 g/sec shows higher maximum pressuriza- 
tion rate and pressure peak than expected from the trend de- 
duced from tests at higher m^g (see Figs. 40 and 44) . This 
is an example of dynamic-burning effect caused by long pre- 
heating of the propellant. 

Calculated pressure delay time (time from igniter onset 
to the instant at which 10% of the maximum pressure is attain- 
ed at the fore-end, t^p) and time to attain maximum fore-end 
pressure are plotted vs Ap/Aj in Fig. 46 for two igniter mass 
flow rates. Both times decrease appreciably when increasing 
m^g (see also Tables 2 and 3). For A p /A t = 1.2 increase of 
mj_g from 11.5 to 19.5 (by 7 0%) decreases i^p from 77.5 to 
37.0 msec (by 52%) . 

2. Effect of Igniter Gas Temperature 

The effective igniter gas temperature (T^g) has a 
strong effect on the ignition delay time and the initial flame 
propagation rate (or the total flame spreading time) and thus 
affects the starting transient in terms of time. Increase of 
T ig by 10% at the datum case (see Table 3) shortens the ig- 
nition time by 24% and the flame spreading time by 14%, and 
has no effect on maximum pressure or maximum pressurization 
rate, because m^g becomes a small fraction of the total 
mass flow rate in the motor in the early stages of the flame 
spreading process. 

3. Effect of the Heat-Transfer Coefficient 


The convective heat transfer correlation affects 
strongly the induction period and flame spreading process. 

An increase of the coefficient a^ in Eq. (III-31) from 
1.56 x 10“ 3 to 1.704 x 10~ 3 (by 9.2%) at the Datum Case 
shortens the ignition delay time by 13% and the flame spread- 
ing time by 9%, and increases slightly the maximum pressur- 
ization rates. If the effect of Ap/A t on the heat-transfer 
coefficient (see Section D, Chapter III) is not considered, 
there is slight reduction in ignition delay and flame spread- 
ing times when Ap/At is increased from 1.2 to 1.5, using the 
same expression for the heat-transfer' coefficient and the 
same igniter gas temperature (see Table 3) . 

4. The Effect of Friction 


The gas dynamic effect alone of friction is to re- 
duce ignition delay and flame spreading times when the fric- 
tion coefficient is increased, because increased friction 
causes higher pressure distribution in the port, and thus 
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increased heat transfer to the propellant. For Ap/A^ = 1.5 
and m^ = 11.5 g/sec, use of 10 times smaller friction co- 
efficient [by multiplying Eq. (III-33) by 0.1] prolongs the 
ignition delay and flame spreading times by about 5%. This 
proves that wall-friction terms have to be included ..n the 
analysis if precise predict ion of the aforementioned times 
is needed. However, the effect of relative surface roughness 
is small and decreases further with increase of the port 
diameter. For example, for the datum case (Ap/A^ = 1.2, 
rtij- = 15.0 g/sec), an increase of the relative equivalent 
sand roughness (e s /du) from 10“3 to 10 -2 shortens the 
ignition delay time by approximately 5% and the flame spread- 
ing tii..e by approximately 2% . 

5. The Effect of Uncertainty in Propellant Properties 

Uncertainty in the propellant thermal properties 
affects the induction and flame spreading periods. Figure 
47 shows the effect of uncertainty of 20% in the propellant 
thermal conductivity (and thermal diffusivity, accordingly) 
on the fore-end pressure transient for Ap/A^ = 1.2 and 
n,ig = 15.0 g/sec. The effect is also seen in Eq. (III-50) . 
The ignition delay and flame spreading times are shorter by 
20% and 16%, respectively, if Xp r is smaller by 20%. The 
propellant ignition temperature (Tp Sf i g ) has, apparently, a 
similar effect. 

The burning rate law affects, obviously, mostly the 
later part of flame spreading process and the entire chamber 
filling period. Therefore, for the same configuration, max- 
imum pressurization rates and pressure levels are determined 
by the function r (Tp^ ,p, pu) . For example, changing the 
erosive burning exponent (6) from 105 to 115 for the datum 
case results in about 5% decrease in the maximum calculated 
fore-and pressure and pressurization rates. 

For the configuration and range of investigation consid- 
ered, the parametric 3tudy may be summarized as follows: 

The ignition delay and flame spreading times decrease 

with: 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


increasing the igniter mass flow rate, m^g,* 
increasing the igniter gas temperature, ; 

increasing A /A. (in part of the range) ; 

r *■’ 

increasing the constant coefficient in the correla- 
tion for heat-transfer coefficient; 

increasing the relative propellant surface roughness; 


FSW" 


* 


' f *„ . 
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( 6 ) 


decreasing the thermal conductivity of the propel- 
lant, Xp r , and 


(7) decreasing the ignition temperature of the propel- 
lant, T . . 

ps, ig 

The maximum chamber pressure and pressurization rate 
increase with: 


(1) increasing the igniter mass flow rate, m ig ; 

(2) increasing the burning surface-to-throat area 
ratio, Aj 3 /A t ; 

(3) increasing the burning rate with respect to a 
reference value at given T . , p, and pu; 

(4) increasing A /A. , and 

p x- 

(5) decreasing the distance from the motor fore-end 
part. 


The dependence of the starting pressure vs time curve 
shape upon various important parameters is shown schematically 
in Fig. 48. 
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CHAPTER VII 

SUGGESTIONS FOR FUTURE STUDIES 

A. Extension of the Analysis 

The next step in the investigation of the starting tran- 
sient of HVT motors is extension of the model and analysis 
developed in this study to various classes of propellants 
and configurations. Certain applications (e.g., use of high- 
ly aluminized propellants) may require modification of the 
analysis and relaxation of some assumption?. However, the 
model is a detailed and fundamental one and any modifications, 
such as using different burning rate laws, inclusion of radi- 
ative heat transfer and consideration of various burning 
surface vs burnt distance functions can be easily introduced. 
Accompaniment of the analysis extension by experimental 
studies or results will verify the validity of the model and 
guide the necessary modifications. Increased versatility of 
the analysis will be very helpful to designers of high-per- 
formance rocket motors. In particular, the techniques should 
be extended to include the radial flow effects in segmented 
motors . 

B. Study of Extreme Cases 

The analysis of extreme cases, such as marginal igniter 
mass flow rate leading to hangfires or misfires, and over- 
sized igniter resulting in high pressure peak, is very impor- 
tant to the proper design and performance of solid-propellant 
rocket motors. Study of the effect of marginal igniter on 
the transient may require consideration of exothermic pre- 
ignition reactions. 5 Analysis of vigorous ignition by over- 
sized igniter, on the other hand, may require consideration 
of dynamic burning contribution. 

The study of both cases as well as the broadening of the 
parametric study presented here are recommended as leading to 
further understanding and ability to control the starting 
transient of HVT motors. 

C. Heat-Transfer Correlation 


\ 



Due to the importance of the heat transfer correlation 
to proper transient analysis and prediction, it is recommended 
to conduct studies for determination of this correlation for 
HVT motors under different practical conditions. The experi- 
mental motor used in this study is an excellent tool for such 
investigations. Experimental heat-flux ifleasuremeDt and perti- 
nent numerical computation techniques 61 ' 68 ' 211-213 are contin- 
uously developing. Detailed knowledge of the fluid dynamic 
flow characteristics is also important and m be obtained by 
optical recording methods. 50 ' 53 Proper evn aation of the 
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effects and phenomena discussed in Section D, Chapter III is 
necessary for determination of the heat transfer to the pro- 
pellant surface. The more is known about the heat transfer 
process, the better the prediction and control of the tran- 
sient become. 


Burning Rate Laws 


The analysis developed in this study considers a burning 
rate which is a function of the initial propellant tempera- 
ture (T • ) , pressure [p(t,x)], heat transfer coefficient 
[h c (t,x7J and mass velocity [p (t ,x) u (t ,x) ] , specific for the 
propellant investigated. One of the conclusions (see Chapter 
VIII) is that improvement of the analysis and prediction as 
applied to HVT motors may be obtained by consideration of 
unsteady (dynamic) burning, as discussed in Section F, Chap- 
ter III and in Appendix C. This calls for a .jroad theoreti- 
cal as well as experimental study of the coupled effect of 
erosive and dynamic burning. Analytically, a burning rate 
law of the following general form is needed: 


r(t,x) = r g [Tpi#p (t,x) ] p p(t,x)u(t,x) ,|£(t,x) ,T pr (t,x,y)J 

Experimentally, the microwave doppier 3hift technique devel 
oped by Shelton^l 1 * and Strand, et al^-> is very promising. 
The direct photographic approach described in Section D, 
Chapter V may be considered too. 
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CHAPTER VIII 
SU MMAR V AND CONCLUSIONS 

A succejsful analytical model was developed to describe 
and analyze the entire starting pressure transient of solid- 
propellant rocket motors with low port-to-throat area ratios, 
resulting in high internal gas velocities. The model was 
formulated in a general way, which enables an easy extension 
of application over a wide range of propellants and motor 
configurations. The flame spreading is coupled to the cham- 
ber gas dynamics through the solid phase heat equation and 
the ignition criterion. Velocity, pressure, and temperature 
variations can now be calculated as a function of time and 
axial distance. 

h suitable numerical method was developed for an effi- 
cient solution of the complicated mathematical model. The 
pre ure transients, calculated by the analytical model are 
in close agreement with the p(t,x) measurements in an exper- 
imental ...otor, utilizing a head -end, pyrogen-type ignition 
system. 

In comparison with previous studies, as reviewed by 
Most and Summerfield ir. Ref. 5 and in thiu work, the essential 
new elements in this study are: 

(1) ability to consider significant pressure, velocity 
and temperature spatial gradients, encountered in 
HVT (High Velocity Transient) motors and their 
variation with time during the starting transient; 

(2) calculation of the propellant surface heating-lo- 
ignition coupled to both temporal and spatial chang- 
es of flow parameters in the chamber during the 
induction and flame spreading periods; 

(3) techniques to account for the strong contribution 
of erosive burning, coupled to the chamber gas 
dynamics ; 

(4) account for axial variation of the port area during 
all three phases of the transient, and 

(5) ability for mo.e precise transient thrust calcula- 
tion and control. 

For port-to-throat area ratios of 1.2, 1.5 and 2.0 an 
agreement within 10% between calculated and measured ignition 
delay times, flame spreading rates, pressure peaks, and times 
to attain the pressure peaks was obtained. For A p /At =1.06 
deviations of up to 15% in time were caused by tho high gas 
velocities and low induction pressures. Better agreement in 
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this case may be obtained by appropriate ignition criterion. 
Maximum pressurization rates as low as 3,000 atm/sec for 
port-to-throat area ratio (A p /A P ) of 1.06 and igniter mass 
flow rate (m^g) of 15 g/sec, ana as high as 7,000 atm/sec 
for Ap/A t = 2.0 and m^g = 19.5 g/sec were measured. The 
highest pressurization rates in a test were obtained in the 
middle of the motor, with slight decrease upstream and large 
decrease downstream of mid-motor. 

Parametric studies have shown that ignition delay and 
flame spreading times decrease with increasing: (1) the ig- 

niter mass flow rate and temperature; (2) A p /A^ between 
1.06 and 1.5; (3) the leading constant in tne correlation 

for heat-transfer coefficient, and (4) the relative propel- 
lant surface roughness. The aforementioned times decrease 
with decreasing the thermal conductivity and ignition tem- 
perature of the propellant. The maximum chamber pressure 
4 nd pressurization rate increase with increasing: (1) the 

igniter mass flow rate; (2) the burning surface-to-throat 
area ratio; (3) the burning rate, (especially the sensi- 
tivity to erosive burning) , and (4) A p /A t , and with decreas- 
ing the d: 3tance from the motor fore-end part. 

When more precise calculation of the starting transient 
is required in particular investigations and applications, it 
can be obtained by using more comprehensive correlations for 
the convective heat-transfer coefficient and the burning rate, 
by taking into account effects, such as increased heat trans- 
fer to rough propellant surfaces, attenuation of heat transfer 
by the "blowing" at the propellant surface due to pre-ignition 
gasification, unsteady burning, resulting from preheating dur- 
ing the induction interval and high pressurization rates dur- 
ing the chamber pressure rise, and turbulent boundary layer 
structure. For very high gas velocities, an ignition crite- 
rion different from that of attaining a critical surface tem- 
perature for ignition is needed. Study of coupled erosive and 
dynamic burning may improve the applicability of the analysis 
presented . 

The experimental system can be readily used to conduct 
various diagnostic tests, pertinent to modern rocket motors, 
e.g., heat-flux measurements, start-stop control, propagation 
of pressure perturbations, and fast liquid-quench extinguish- 
ment . 


The results of this study are useful to designers of 
high-performance solid-propellant motors. In particular, by 
using the techniques developed in this work and the results 
obtained, the designer can prescribe the upper limit of volu- 
metric loading density, that can be used, -nd predict more 
accurately the pressure-thrust-time motor performance during 
the starting transient. 
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Table 1 

Baseline Values 

Propellant properties : 

Composition: 20% PBAA-EPON 828 

80% AP (30% 15 micron & 70% 180 micron) 

r 0@p=68 atm = °* 804 cm/sec = °' 32 in/sec 
n = 0.4 

_3 

Xp r = 0.9 x 10 cal/cm-sec-°K 

3 

Ppr =1.6 g/cm 
a = 0.1875 x 10" 2 cm 2 /3ec 

p3T 

W = 22.01 g/g-mole 
y = 1.24 

Cp = 0.4665 cal/g-°K 

R = 3852 gf-cm/g-°K 

T . = 298°K 
Pi 

T ■ = 7 00°K 

ps r lg 

T f = 2225°K 


k = 5.72 am -°K/cal 
8 = 105 


e„ = 0.001 cm 
s 


Motor parameters : 

x_ = 6.2 cm 
P 

b = 5.08 cm 


A . = 1.61 cm‘ 

Pf i 

L* = 62.7 (A p /A t ) cm 


x_, = 62.7 cm 

E 

d. . = 1.02 cm 
h,i 

P . = 6 . 3 5 cm 
w,i 
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Table 4 

Typical Mach and Reynolds Numbers (based on diameter) during 
the different stages of the transient for the four port-to- 

throat area ratios tested and m. = 11.5 g/sec 

J-9 


Vt 


1.06 



1.5 


2.0 

It 

II 


Period 

Fore-end 
Mach No. 

Aft-erd 
Mach No. 

Fore-end 
Reynolds No 

Aft-end 
. Reynolds 

induction 

0.58 

0.42 

11,500 

17,000 

flame-spread ing 

0.080 

0.76 

10,000 

100,000 

chamber filling 

0.016 

0.76 

10,000 

320,000 

peak pressure 

0.017 

0.76 

10,900 

350,000 

induction 

0.55 

0.39 

11,500 

17,500 

flame spreading 

0.090 

0.60 

10,000 

100,000 

chamber filling 

0.015 

0.60 

10,000 

330,000 

peak pressure 

0.016 

0.60 

11,000 

360,000 

induct ion 

0.52 

0.35 

11,700 

19,000 

flame spreading 

0.090 

0.44 

10,000 

80,000 

chamber filling 

0.010 

0.44 

10,000 

360,000 

peak pressure 

0.011 

0.44 

11,100 

400,000 

induction 

0.45 

0.28 

12,000 

21,000 

flame spreading 

0.060 

0.31 

10,000 

50,000 

chamber filling 

0.009 

0.31 

12,000 

400,000 
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ig. 3 Cross-section of experimental motor 
showing port configuration. 
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LOCATION CF MAXIMUM 
HEAT TRANSFER 



i (a) duct of relatively small diameter 
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Pig* 5 Jet expansion from a sonic convergent nozzle into 
a duct (from Ref. 53) . 
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Fig. 6 Non-erosive burning rate law of propellant PBAA COBM at initial 
temperature of 298°K. 
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Exploded view of the experimental window motor 
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Overall view of the assembled gaseous ignition system mounted on a test stand 






Typical variation with time of the pressure at motor fore-end and 
igniter mass flow rate, showing the instant of igniter nozzle unchoking 
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Propellant cast in trays and on Plexiglas window for photographic observation 
of ignition and receding propellant surface during burnina. 
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Fig. 22 Experimental starting pressure transients recorded during a test 
with A /A. = 1.5. 
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TIME FRON ONSET OF IGNITER, t, MSEC 
Experimental starting pressure transients recorded during a test 





PRESSURE AT STATION 1 (FORE-END) , p, , KGF/CM 


L2" 



TIME FROM ONSET OF IGNITER, t, MSEC 

Fig. 25a Comparison of measured and calculated pressure 

vs time traces at station 1 (fore-end) for 

A /A. = 1.5. 
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PRESSURE AT STATION 3 (MID-MOTOR), p_, KGF/CM 


RUN D-9 
A /A t =1.5 



Fig. 25b Comparison of measured and calculated pressure 
vs time traces at station 3 (mid-motor) for 
A / A . = 1.5. 









PRESSURE AT STATION 5 (AFT-END) 


RUN D-9 


V A t = 1 - 5 

m. =17.9 G/SEC 
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T. = 2340°K 
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TIME FROM ONSET OF IGNITER, t, MSEC 

Comparison of measured and calculated pressure 
vs time traces at station 5 (aft-end) for 
A p /A fc » 1.5. 
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Fig. 26a Comparison of measured and calculated pressure 
vs time traces at station 1 (fore-end) for 
A p /A fc -1.2. 



PRESSURE AT STATION 3 (MID-MOTOR), p, , KGF/CM 



TIME FROM ONSET OF I (NITER, t, MSEC 

Fig. 26b Comparison of measured and calculated pressure 
vs time traces at station 3 (mid-motor) for 

V A t ■ 1 
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30 40 50 60 70 80 

TIME FROM ONSET OF IGNITER, t, MSEC 

Fig. 26c Comparison of measured and calculated pressure 
vs time traces at station 5 (aft-end) for 
A /A t -1.2. 






PRESSURE AT STATION 1 (FORE-END) , p. , KGF/CM 
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EROSIVE BURNING 

CALCULATED COMPLETION 
OF FLAME SPREADING 


TIME FROM ONSET OF IGNITER, t, MSEC 

Fig. 27a Comparison of measured and calculated 
pressure vs time traces at station 1 

(fore-end) for A /A. = 1.06. 
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PRESSURE AT STATION 3 (MID-MOTOR), p , KGF/CM 






PRESSURE AT STATION 5 (AFT-END) , p_ , EGF/CM 



30 40 50 60 70 80 


TIME FROM ONSET OF IGNITER, t, MSEC 

Fig. 27c Comparison of measured and calculated pressure 
vs time traces at station 5 (aft-end) for 
A p /A t = 1.06. 
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HEAT FLUX AT THE PROPELLANT SURFACE , q, CAL/CM -SEC 
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TIME FROM ONSET OF IGNITER, t, MSEC 


Fig. 35 Calculated heat flux at the propellant surface 
vs time at three different positions along the 
motor . 



PROPELLANT SURFACE TEMPERATURE 
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Fig. 36 Calculated propellent surface temperature vs 
time at three different positions along the 
port. The calculated surface heat flux is 
shown in Fig. 35. 
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LOCATION OF THE FLAME FRONT, x - x p , CA 

Fig. 38 Measured and calculated ignition flame 
front locations vc time, with flame 
spreading velocities noted. 
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INITIAL FLAME SPREADING RATE 







MEASURED MAXIMUM PRESSURIZATION RATE AT MOTOR FORE -END, KGF/CM -SEC (x 10 



IGNITER MASS FLOW RATE, m ig , G/SEC 

Fig. 40 Measured maximum pressurization rates at motor 

fore-end (station 1) vs igniter mass flow rate 

for different A /A. . 
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TIME PROM ONSET OF IGNITER, t, MSEC 

Calculated transient variation of the fore-end and 
aft-end static pressures, aft-epd stagnation pressure 
and ratio of aft-end stagnation to fore-end static 
pressure. 


Fig 
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PORT-TO-THROAT AREA RATIO (A p /A fc ) 

Fig. 43 Calculated ratios of maximum static and stagna- 
tion pressures at motor aft-end to maximum 
motor fore-end pressure vs port-to-throat area 
ratio. 
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Fig. 47 Effect of uncertainty in propellant thermal conductivity on the 




Fig. 48 Schematic representation of the effect of various parameters 
on the shape of the starting pressure transient of HVT motors 
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APPENDIX A 

NONDIMENSIONALIZATION AND ORDER OF MAGNITUDE 
ANALYSIS OF THE CONSERVATION EQUATIONS 
1. Main Reference Quantities 


Nondimensionalization of the conservation equations is 
performed by normalizing all variable parameters with respect 
to well-defined reference quantities. The main reference 
quantities used for the nondimensionalization of the conser- 
vation equations [Eqs. (IIJ-1) to (III-3)] are the following: 


(1) Reference time, t*, defined as the average resi- 
dence time of the gas in the motor chamber; 

t * _ instantaneous mass of gas in the chamber 

mass rate of gas ejected out of the chamber 

Assuming that the representative static pressure 
and temperature in the motor can be approximated 
by the stagnation pressure and temperature at the 
nozzle entrance, (p as and T as , respectively) , 
the following expression for t* may be obtained: 


t* 


L 

r^c* 



(A-l) 


where: L is the effective length of the chamber 

(the distance between the fore-end of 
the propellant section and the nozzle 
entrance) 

c* is the characteristic velocity 
and T is a function of y, defined in the 
Nomenclature. 


(2) Reference length, x*, defined as the effective 
length of the motor chamber (x* = L) ; 

(3) Reference gas velocity, u*, defined as the aver- 
age gas velocity in the chamber port 


u* 



r 




(A-2) 


(4) Reference pressure, p*, defined as the maximum 
pressure obtained during the transient period at 
the motor fore-end (p* « p.^ ) ; 
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(5) Reference gas temperature, T*, defined as the 
adiabatic flame temperature of the solid pro- 
pellant (T* = T f ) ; 

(6) Reference density, p*, determined by the refer- 
ence Pressure and temperature, and 

(7) Reference burning rate, r*, defined as the non- 
erosive steady-state burning rate at the refer- 
ence pressure (r* « ap* n ) 

All configuration (i.e., geometrical) parameters are 
nondimensionalized with respect to their initial values. 


All reference quantities are denoted by the superscript 
and all nondimensionalized variables are denoted by the 
sign " A ". According to this convention 


£ “ ; u ■ -T ; T * ? £ 


t = pf ; & = j-f •' etc. 



Derivatives with respect to 

3p _ p* 9b 3u u* 
E-Faf f It 


time become 
etc . 


(A-3) 


(A-4) 


Derivatives with respect to distance become 


3p _ p* 36 . 3u _ u* 3d 
« ‘ x 7 3$ ' 3x ” x 7 31F 


etc . 


(A— 5 ) 


Second derivatives with respect to distance become 
3 2 T T* 3 2 T 3 2 u _ u* 3 2 fl 

K ~ K K ; K - K 5 " 

3x z x*^ 3iT 3x^ x** 3iT 


2. Nondimensionalization of the Continr .ty Equation 

The continui y equation, Eq. (111—1 ) , is nondimension- 
alized according to Eqs. (A-2) to (A- 5) to become 


P*U* 

x^“ 



+ 



p Sr r * b * fB 


pr 



(A-7) 
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Dividing Eq. (A-6) by p*u*/x* the nondimensional- 
ized continuity equation is obtained: 



(A-8) 

where an order of magnitude characteristic of the physical 
model is denoted under the terms. 


The last term on the right-hand side of Eq. (A-8) is 
10 2 times smaller than the other terms for the usual oper- 
ating pressures, flame temperatures, gas molecular weights, 
and propellant densities encountered in solid-propellant 
rocket motors (20 to 80 atm; 2000 to 3000°K; 20 to 30 g/g- 
mole, and 1.5 to 2.0 g/cm , respectively). It is therefore 
neglected in the analysis. The neglected term expresses 
the rate of mass accumulation in the free volume created 
by the propellant surface regression. After this simpli- 
fication, the resulting dimensional continuity equation is 
presented by Eq. (III-4) . 


3. Nondimensionalization of the Momentum Equation 


Substituting Eq. (III-l) into the momentum equation, 
Eq. (III-2) , the latter becomes (using the expression for 
normal stress) 

fP. 


du 


p Tt + pu 37 + 


2u , 3p p pr rl5 , 4 
— + g-*f“ ■ - — u + jy 


ifu 1_ 3u 

17 A p 8 * 8 * 


w 2 

3a^ u 

(A-9) 

where y may be considered to include the eddy viscosity. 
Nondimensionalization of Eq. (A-9) yields 


p* (u*) 


*\ 2 


f.aal . p*<u*) 2 L„3<ti . YcrR*T*_.ae 
+ ( P< %J + YX* ° ~5t 


P*r*b* 6 

- 

p p 


, 4 #n u* |a 2 a . 1 ^P 3fl| f * P w * . 2 

+ 5^ S3T 33 ■ “Sa* 5 {u > (a ‘ 10) 

1 *> 2 . 


Dividing throughout by p*(u*) /x* and denoting 


Re 3 


p * u *d* 
— 7 * — 


and (M*)» (u*) 2 /ygR*T* 


--r 3S 


«?• - _ 
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I 




10~ 4 Iff 2 11 1 10“ 2 10 1 

' (A-ll) 

where an order of magnitude characteristic of the physical 


i 

I 


model is denoted under the terms. 

The second term on the right-hand aide of Eq. (A-ll) , 
expressing the viscous forces between gas molecules is 
several orders of magnitude smaller than other terms. 
Therefore it is neglected in the analysis. The torm des- 
cribing the wall friction forces acting on the fluid is 
one order of magnitude smaller than the other terms and is, 
therefore, kept in the analysis. Other reasons for keeping 
the wall-friction term are described in Chapter III , Sec- 
t : on C. After the aforementioned simplification, the 
resulting dimensional momentum equation is presented by 
Eq. (III-5) . 


Nondimensionalization ol the Energy Equation 



Substituting the original continuity and momentum 
equations (Eqs. (III-l) and (III-2)] into the energy equa 
tion, Eq. (III-3) , the latter becomes (in terms of tbs 
static gas temperature) 
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Nondimensionalization of Eq. (A-12) yields 


p*c*T*u* 

— 1 2=*— P6. 


w + 


3T' 


Y-l P * C J T * U * 
Y x* 


$*♦•3 


T* ,3 2 t P^ r r*b*c;T* f fM 

' x 1^*5? = *T 


- T + g- 


4q £ . A*T* X 3t 3A p , 4 M*(u*) 2 . 
3 ££ 

, 2 f*p*{u*r 1 w* A 3 

+ T* 

p 


2 r 9a) 2 


air 


(A-13) 


Assuming the heat transfer to the port walls to be 
convective and using the Dittus-Boelter correlation for the 
heat-transfer coefficient (Nu. = a, Re 0 • 8pr0 *4) , the former 
may be expressed by the following equation 

<J t - Hu a a^ (T - V 


“ Nu S ar<* - V 

n h 


(A-14) 


Substituting Eq. (A-14) into Eq. (a- 13) and dividirg 
through by p*c*T*u*/x* yields the following form of the 
nondimensional energy equation and the corresponding order 
of magnitude of the different terms: 




1 


1 

Re*Pr* 




b*x* 



fi f Be 

P F. r ^ . £ 






1 


10 2 ltf' 10 2 


111 1 


oio 
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lO -1 10~ 5 ID -2 1 10" 2 lCf 1 10 2 1 

(A-15) 

where St* is the reference Stanton Number (St* = Nu*/ 
Re*Pr*) . 

The following terms in Eq. (A-15) are several orders 
of magnitude smaller than the other terms: (1) the third 

term on the left-hand side and the third term on the right- 
hand side, which describe the axial heat conduction between 
gas molecules; (2) the fourth term on the right-hand side, 
which describes the viscous dissipation. Therefore, these 
terms are neglected in the analysis. The term describing 
the rate of work performed on the fluid by the wall fric- 
tion forces is one order of magnitude smaller than the 
other terms and is kept in the analysis. After this sim- 
plification the resulting dimensional energy equation is 
presented by Eq. (III-6) . 
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APPENDIX B 

THE AXIAL VARIATION OF PRESSURE , VELOCITY AND 
TEMPERATURE IN THE MOTOR PORT 

The effects of the various processes taking place dur- 
ing the pressure transient on the axial variation of pres- 
sure, velocity and temoerature are shown in the following 
paragraphs . 


Equations (III-7) and (III— 9) can be solved for the 
par'-ial space derivatives of pressure and gas velocity to 
yield the following expressions: 

r 

3p _ M Z J 1 3p _ yp 9u 

3x , ..2 1 u at 2 at 


3A 

+ IE _J 

A 3x 
P 


VARIATIONS W.R.T. TIME 
(UNSTEADINESS EFFECTS) 
P„rb _ f T , 


- yR- 


A u 
P 


T , , f . Y-1 m 2 

u x + T" + T 41 


EFFECT OF AXIAL 
PORT AREA VARIATION 


EFFECT OF PROPELLANT GASES 
ENTERING THE STREAM 


(y-i)j p 

, 1 c w „ 

A u •*£. 
P 


- X — p F ( Y “ 1 ) M 2 + 1 ► 

P L J 


(B-l ) 


EFFECT OF HEAT 
TRANSFER TO PORT WALLS 


EFFECT OF FRICTION 
LOSSES 


JL M2 au 1 3£ 

2 u 3t YP St 


A p 5x 


RP P frb * U±±M 2 + J* 

Ap Pi 2 T 


p w <?* 
Y Ap P 


w 

: — fM u 


(B-2) 
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From Eq. (III-8) an expression for the partial space 
derivative of gas temperature is obtained: 

a a 

3T _ 1 9T _ , ,,T 3u _ (y-1)T 

3x u 9t u 9x A 5x 

P 


+ yR- 


p pr rb 

P u 




- (Y-l)J 


w T 


c A pu'H 
P 


+ Yty-P P w fu „2 

2 A_ 


(B-3) 


The effects of the various parameters on the axial 
variation of pressure, temperature and velocity are clear- 
ly expressed in Eqs. (B-l) througn (B-3). The pressure 
variation along the port [Eq. (B-l)] is of particular 
interest. For the subsonic flow existing in the motor, 
increases of cross-sectional port area and heat loss to 
the port walls tend to increase the pressure with distance 
from the fore-end, whereas increases of added mass and 
friction cause axial decrease of pressure. During the 
induction period the unsteadiness effects are negligible 
as compared with the effects of other terms, except for a 
very short time (1 to 4 msec) at the igniter onset. The 
pressure distribution (and the distribution of temperature 
and velocity as well) is determined by the heat flux to 
the port walls, the wall friction, and the initial port 
area distribution. During the flame spreading period and 
the early part of chamber filling interval, however, the 
flow is largely unsteady and the effects of the terms con- 
taining time derivatives [first two terms on the right-hand 
side of Eq. (B-l)] become significant. The dominant term 
is that expressing the effect of added mass [third term 
on the right-hand side of Eq. (B-l)]. Typically, for 
Ap/A t = 1.2 the combined magnitude of unsteadiness terms 
for flow over the unignited part of the port during flame 
spreading is between 10 and 20% of that of the heat-trans- 
fer term. Immediately following the completion of flame 
spreading the magnitude of unsteadiness terms at the fore- 
end of the propellant section is about 10% of that of the 
mass-addition term. Thereafter the effect of unsteadiness 
decreases and that of mass addition increases. 
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APPENDIX C 

COMPARISON OF SOLID PHASE RELAXATION TIMES WITH 
THE CHARACTERISTIC TIME OF MOTOR PRESSURIZATION 


The phenomenon of dynamic burning following ignition of 
the solid propellant due to preheating is discussed in Sec- 
tion F, Chapter III. An estimate of its importance in view 
of the characteristic time of pressure change during the 
transient (t ) is given in the following paragraphs. 


The characteristic thickness of preheated layer at the 
time of ignition (6 . ) , assuming constant heat flux to an 
inert propellant suPrace can be estimated from the expres- 
sion 


ph 


T — T X(T — T) 

ps , ig pi _ pr ps,ig pi 

OT_/ay)„ = J 


(C-l) 


pr' 


y=0 


The time to burn through the preheated layer (t . ) is 
given by p 


T ph ^ph^ r 


(C-2) 


where r is an average burning rate evaluated at average 
pressure between the induction-period pressure and the maxi- 
mum pressure. 

Consider the following two definitions of characteristic 
time for motor pressurization: (1) the total time of pressur- 

ization, x . , from the instant of pressure rise following 
the induct?6n-interval level to the attainment of maximum 
pressure; and (2) the pressurization time interval corres- 
ponding to maximum pressure rise, x . The latter is given 
by the relation p ' 


p max p ind 

T p,m = (9P/3t) max 


(C-3) 


where P ind is the pressure level during the induction inter- 
val. 

For a typical test with the experimental motor at A /A. 
=1.2, x .=40 msec and x = 10 msec (at the fore end or 
the propillant section) . p ' 

For the baseline values listed in Table 1 and an average 
constant value of 40 cal/cm -sec for q. (see Fig. 35) ex- 
pression (C-l) yields a value of 0.009 cm for 6 . * if 
burning rate in Eq. (C-2) is evaluated as the non-erosive 
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2 

burning rate at an averaqe pressure of 20 kgf/cm (see Fig. 

6), then the time to burn through the preheated layer, 
is 18 msec. An estimation of the relaxation 2 time for 
solid phase steady-state burning (t = u /r ) for the same 
burning rate yields a value of 7.8 fisec.^ 

Comparing these values of t . and x with the afore- 
mentioned typical values for x and x it may be con- 

cluded that, as far as the entir4 m pressuri?&tion time is con- 
cerned, the effect of dynamic burning due to preheating is 
small. For the limited time of fast pressure change 

the effect is significant. However, for properly designed 
motors the period of large dynamic burning effects is over 
before the pressure overshoot occurs . ^6 • ^16 Thus the dyn- 
amic burning is not expected to significantly increase P max * 
As mentioned in Chapter VII, combined consideration of ero- 
sive and dynamic burning may improve the analysis and extend 
its application. 
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APPENDIX F 

LISTS OF COMPUTER PROGRAM SUBROUTINES 

AND INPUT DATA CARDS 

1. Description of Subroutines 

The computer program consists of a main program and 
thirteen subroutines. 

Subroutines LBCUOP and RBC calculate the left- 
boundary and right-boundary values of the flow parameters, 
respectively. 

Subroutine TSCAL calculates the propellant surface 
temperature by integrating Eq. (III-50) . 

Subroutine TPSABL calculates the propellant surface 
temperature during a decomposition-gasification period 
(when assumed) by integrating an equation similar to Eq. 
(III-50) obtained by the integral method. 

Subroutines TRID, DECOMP, SOLV, GSOLV, ASOLV, VSOLV, 
REPL2 , REPLl , and SING solve the block-tridiagonal matrix 
of the finite-difference equations by decomposing it into 
two bidiagonal matrices. 

2 . Input Data Cards 

There are nine input data cards, placed at the end of 
the computer program deck. Each card contains up to seven 
input parameters, placed in consecutive groups of ten 
columns. The format of each parameter is F10. Therefore 
the number may be specified either as X.XXXXE+OX or 
XXXX.XXXX, as required. The designation in the following 
list is the notation used in the computer program. 

The First Card contains the following input information: 

1) XP = position at the entrance to the propel- 

lant section, cm (x in the text); 

Er 

2) XG = position at the aft-end of the propel- 

lant slab, cm; 

3) XE = position at the entrance to the motor 

nozzle, cm (x £ in the text) ; 

2 

4) API = initial port area, cm ; 

5) API AT = initial port-to-throat area ratio 

(A /A. in the text) ; 

P 
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6) BP = burning perimeter, cm (b in the text); 

7) WPI - initial wetted perimeter, cm (P 

in the text) . w ' 1 

The Second Card contains the following input informa- 
tion : 


1) W = molecular weight of the combustion gases, 

g/g-mole (W in the text) ; 

2) GAMA = ratio of specific heats (y in the 

text) ; 

3) ROPR = density of the solid propellant, 

g/cm 3 (p pr in the text) ; 

4) FKPR = thermal conductivity of the solid 

propellant, cal/cm-sec-°K (A in 
the text) ; P r 

5) ALPHAP = thermal diffusivity of the solid 

propellant, cm2/sec (a in the text) 
(c is calculated in ?fie program 
froA inputs 3, 4, and 5) 

6) Leave blank; 


7) ROUGH = equivalent relative roughness of the 

port walls (e /d. in the text) . 
s n 

The Third Card contains the following input informa- 
tion: 


1) TF = adiabatic flame temperature of the solid 

propellant, °K (T^ in the text); 

2) TPSCRI = critical propellant surface tempera- 

ture for start of decomposition- 
gasification, °K ; 

3) TPSBRN = propellant surface temperature, at 

which propellant ignition occurs, 

(T . in the text); 

ps,ig 

4) A = pre-exponential factor in the non-erosive 
burning rate law, ap*i (a in the text) ; 

BREXP = pressure exponent in the non-erosive 
burning rate law (n in the text) ; 


5) 
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3 

6) EBC = erosive burning constant, cm -°K/cal 

(k in the text) ; 

7) EBEX = erosive burning exponent (3 in the 

text) . 

The Fourth Card contains the following input informa- 
tion: 

1) TCOR = time of the first corner in the multi- 

linearized m. vs t plot (see Fig. 14 ) , 
lg 

sec ; 3 

2) TSS = time, at which steady-state igniter 

operation is reached (second corner in 
Fig. 14 ) , sec; 

3) TCUT = time, at which linear igniter cut-off 

starts (third corner in Fig. 14), sec; 

4) TOFF = time, at which igniter is completely 

turned off, sec; 

5) MIGN(O) = initial igniter mass flow rate 

(at t - 0) to facilitate the 

numerical calculation, g/sec, 

(m. . in the text); 

lg , l 

6) MIGIJC = igniter mass flow rate at time TCOR, 

g/sec ; 

7) MIGNSS = steady-state igniter mass flow rate, 

g/sec. 



Time , sec 
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The Fifth Card contains the following input informa- 
tion: 


1) TGI = initial gas temperature in the motor, 

°K (T^ in the text) ; 

2 

2) PCI = initial pressure in the motor, gf/cm 

(p^ in the text) ; 

o 

3) PAM = ambient pressure, gf/cm (p in the text); 

3 . 

4) TPI = initial propellant temperature, °K 

(Tpi in the text) ; 

5) TIGN * effective mean temperature of the 

igniter gas, °K (T^ g in the text) ; 

6) UCI = initial gas velocity in the motor, 

cm/sec (u^ in the text) . 

The Sixth Card contains the following input informa- 
tion: 

1) NITER * number of iterative predictor calcu- 

lations; 

2) NPRINT = number of time-step calculations, 

after each of which an output is 
printed out; 

3) NPUN = number of time-step calculations, af- 

ter each of which an output is punched 
out for plotting and restart; 

4) NSPAVE = number of time >p calculations, 

after each c: \-r h space averaging 
takes place; 

5) NBL - number of tim e-s-^p -a ■ aticvs, after 

each of which balai . ; che conserva- 

tion and governing . ons is perform- 
ed and printed out 

The Seventh Card contains the following input informa- 
tion: 


1) DELTAX * spacewise step for the numerical 

calculation (must result in a whole 
number of increments), cm (Ax in the 
text) ; 

2) DELTAT ■ time step in the numerical calcu) 

tion, sec (At in the text) ; 
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3) TI = initial time of a computer run, or a 

restart run, sec; 

4) TMAX = maximum time of a computer run, or 

a restart run, sec; 

5) EPSLON = a small number in Eq. (III-51 ), 

°K (e in the text) ; 

6) DDRG = a factor, by which the expression for 

the friction coefficient [Eq. (III-33 )] 
is multiplied (for parametric studies) ; 

7) DDHC = a factor, by which the expression 

for the heat-transfer coefficient 
[Eq. (III-31)] is multiplied (for 
parametric studies) . 

The Eighth Card contains the following input informa- 
tion: 


1) TSAVl = time, before which the space-averaging 

weight SAVMIN (see below; is used, 
sec; 

2) TSAV2 * time, after which the space-averaging 

weight SAVMAX (see below) is used, 
sec ; 

3) SAVMIN = space-average percentage weight of 

the values at the neighboring me^h 
points specified for the time before 
TSAVl ; 

4) SAVflAX * space-average percentage weight of 

the values at the neighboring mesh 
points specified for the time after 
TSAV2 . Between the times TSAVl and 
TSAV2 linear increase n the space- 
average weight is specified. 

The Ninth Card contains the following input informa- 
tion : 


1) TTT1 = time, before which the Crank-Nicolson 

parameter has the value THETAN (see 
below) , sec; 

2) TTT2 = time, after which the Crank-Nicolson 

parameter has the value THETAX (see 
below), ■~'!c; 
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3) THETAN = value of Crank-Nicolson parameter, 

specified for the time before TTTl; 

4) THETAX = value of Crank-Nicolson parameter, 

specified for the time after TTT2 . 
between the times TTTl and TTT2 
linear increase in the parameter is 
specified 


i 




Table FI 




List of Input 

Values for 

Datum 

Case 

XP 

2 

6.16 

TGI 

= 

298.0 

XG 

= 

54.42 

PCI 

= 

1.0331 x 10 

XE 

= 

62.04 

PAM 

= 

1.033 x 10 3 

API 

= 

1.613 

TP I 

= 

298.0 

API AT 

= 

1.20 

TIGN 

= 

2410.0 

BP 

= 

5.08 

UCI 

s 

727.0 

WPI 

= 

6.35 

NITER 

= 

1.0 

W 

= 

22.01 

NPRINT 

= 

20.0 

GAMA 

= 

1.24 

NPUN 

= 

400.0 

ROPR 

= 

1.60 

NS PAVE 

= 

1.0 

FKPR 

= 

0.0009 

NBL 

= 

200.0 

ALPHAS 

= 

0.001875 

DELTAX 

= 

2.54 

ROUGH 

= 

0.001 

DELTAT 

= 

0.000025 

TF 

= 

2225.0 

TI 

= 

0.0 

TPSCRI 

s 

690.0 

TMAX 


0.080 

TPSBRN 

= 

700.0 

EPSLON 

= 

0.10 

A 

= 

0.00927 

DDRG 

= 

1.0 

BREXP 

= 

0.40 

DDHC 

= 

1.0 

EBC 

= 

5.72 

TSAV1 

= 

0.01 

EBEX 

= 

105.0 

TSAV2 

= 

0.02 

TCOR 

= 

0.002 

SAVMIN 

= 

0.06 

TSS 

= 

0.012 

SAVMAX 


0.06 

TOUT 

= 

0.150 

TTT1 

= 

0.01 

TOFF 

= 

0.220 

TTT2 

= 

0.02 

MIGN(O) 

= 

1.055 

THETAN 

= 

0.60 

MIGNC 

= 

12.0 

THETAX 

= 

0.60 

MIGNSS 

= 

15.0 






